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Abstract 

The  use  of  elgenstructure  assignment  techniques  has 
received  wide  attention  as  a  tool  for  designing  flight 
control  systems  for  aircraft  with  multiple  control  surfaces. 
One  drawback  for  using  this  technique  is  a  lack  of  handling 
quality  guidelines  to  apply  when  selecting  the  eigenvalues 
and  eigenvectors  of  the  closed-loop  system.  This  lack  of 
specific  elgenstructure  requirements  means  that  some 
uncertainty  will  remain  as  to  whether  the  augmented  control 
system  will  meet  the  MIL-F-8785C  specifications. 

.  Therefore,  development  of  a  method  for  choosing  the 
desired  elgenstructure  of  the  augmented,  closed-loop  system 
which  would  meet  the  handling  qualities  specifications  was 
examined.  This  method  consisted  of  forming  an  "optimal" 
plant  matrix  which  possessed  desirable  dynamic  characteris¬ 
tics  and  performing  a  spectral  decomposition  of  this  matrix. 
The  resulting  elgenstructure  was  used  as  the  desired 
eigenvalues  and  eigenvectors  during  the  full-state  feedback, 
elgenstructure  assignment  process.  The  resulting  feedback 
gain  matrix  was  used  in  the  control  system. 

As  an  example,  this  process  was  performed  on  a  model  of 
the  X-29A  using  the  canard,  flaperon,  and  strake  flap  control 
surfaces.  The  resulting  augmented  system  was  evaluated  using 
the  Neal-Smith  pilot-model  analysis  and  also  using  an  X-29A 
man-ln-the-loop  simulation.  The  results  show  that  the  method 
is  very  promising,  although  care  must  be  taken  that  all 
anticipated  control  system  dynamics  are  considered  when 
forming  the  optimal  A  matrix. 
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APPLICATION  OF  EIGENSTRUCTURE  ASSIGNMENT  TECHNIQUES 
IN  THE  DESIGN  OF  A  LONGITUDINAL  FLIGHT  CONTROL  SYSTEM 


I .  Introduction 

Background 

Because  of  the  Increasing  demands  being  placed  on 
today's  fighter  aircraft,  each  aircraft  system  must  be 

exploited  to  Its  fullest  potential  in  order  to  yield  superior 

■ 

performance  and  maneuverability.  Engineers  within  the 
aerospace  community  are  continually  searching  for  Innovative 
designs  which  will  prove  superior  to  present  and  future 
threats.  One  of  the  consequences  of  trying  to  Increase  the 
maneuverability  of  our  aircraft  Is  that  Inherent  aerodynamic 
Instabilities  are  being  tolerated  by  aircraft  designers  in 
order  to  attain  the  quicker  maneuvering  response  which  static 
Instabilities  allow. 

An  unstable  aircraft  demands  a  closed-loop  flight 
control  system  to  maintain  steady-state  flight.  Pilot 
reaction  time  and  work-load  restrictions  demand  that  the 
aircraft  exhibit  predictable  flying  characteristics,  at  least 
from  the  pilot's  perspective.  In  addition,  multiple  control 
surfaces  per  aircraft  control  axis  are  being  used  to  decouple 
the  characteristic  modes  or  to  allow  for  controlled  flight  In 
the  post-stall  regime  f 1 : 27, 33 J .  Again,  a  flight  control 
system  is  needed  to  deflect  multiple  control  surfaces  in  the 
correct  proportion  so  that  the  aircraft  responds  as  the  pilot 
has  commanded  through  the  cockpit  controls. 

The  use  of  feedback  to  augment  a  control  system  Is  a 
well-known  and  accepted  practice.  It  is  done  for  a  number  of 
reasons  including:  Improving  the  stability  of  the  system, 
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reducing  sensitivity  to  modeling  inaccuracies,  and  changing 
the  response  behavior  of  the  system.  In  the  design  of  an 
aircraft  flight  control  system  for  an  aerodynamlcally 
unstable  aircraft  for  which  only  wind  tunnel  or  other 
empirical  data  exists,  all  three  of  these  reasons  for 
feedback  are  pertinent. 

The  application  of  eigenstructure  assignment  techniques 
using  full-state  or  output  feedback  has  been  developed  and 
used  in  the  last  ten  years  to  fulfill  the  three  previously 
mentioned  requirements,  especially  with  the  emergence  of 
fly-'by-wire  control  systems.  Eigenstructure  assignment  also 
has  advantages  over  classical  design  techniques  for 
Multlple-lnput/Multlple-Output  (MIMO)  systems.  While 
application  of  the  eigenstructure  assignment  techniques  has 
progressed  rapidly,  there  still  exists  some  difficulty 
relating  the  choice  of  the  desired  eigenvalues  and 
eigenvectors  to  the  requirements  of  aircraft  handling 
qualities.  In  the  final  analysis,  the  complete,  augmented 
closed-loop  system  must  behave  in  a  manner  acceptable  to  the 
pilot . 

The  MIL-F-8785C  document  "Military  Specification-Flying 
Qualities  of  Piloted  Airplanes"  [21,  along  with  its  handbook 
"Background  Information  and  User  Guide  for  MIL-F-8785C"  [31, 
currently  contain  the  guidelines  used  by  aircraft  control 
system  engineers  to  help  insure  that  the  resulting  system 
will  be  acceptable  to  the  pilot.  However,  the  specifications 
of  MIL-F-8785C  were  written  when  the  classical  frequency 
domain  methods  of  control  theory  were  used  to  design  flight 
control  systems.  The  control  engineer  using  modern  control 
theory  cannot  directly  design  the  entire  eigenstructure  by 
these  same  specifications.  One  is,  therefore,  somewhat  at  a 
loss  when  choosing  eigenvalues  and  especially  eigenvectors. 


for  the  augmented  control  system.  In  addition,  the  draft  for 
the  new  flying  qualities  standard  does  not  directly  address 
eigenvalues  or  eigenvectors  [41. 

This  lack  of  specific  elgenstructure  requirements  means 
that  despite  the  fact  that  elgenstructure  assignment  is  a 
very  powerful  tool,  some  uncertainty  will  remain  as  to 
whether  the  augmented  control  system  will  meet  the 
MIL-F-8785C  specifications.  Usually,  this  uncertainty  leads 
to  an  iterative  procedure  Involving  analyzing  a  prospective 
control  system  using  a  manned  simulation  or  pilot-model 
analysis  tool,  and  then  modifying  the  control  system  based  on 
the  results  of  the  developmental  studies.  Ideally,  given  a 
reasonably  accurate  aircraft  model,  one  would  like  to  be  able 
to  make  only  one  iteration  in  the  design  of  a  flight  control 
system.  This  means  having  the  capability  of  picking  a 
realizable  elgenstructure  which  will  yield  satisfactory 
handling  qualities  the  first  time. 

Problem 

Development  of  a  method  for  choosing  the  desired 
elgenstructure  of  the  augmented,  closed-loop  system  which 
would  meet  the  handling  qualities  specifications  and, 
furthermore,  achieving  this  elgenstructure  through  the  use  of 
elgenstructure  assignment  would  eliminate  the  uncertainty  of 
using  this  design  technique.  This  entails  examining 
MIL-F-8785C  and  the  associated  Handbook  and  identifying  key 
parameters  which  can  be,  either  directly  or  indirectly, 
translated  into  a  specification  of  the  elgenstructure.  This 
method  should  be  proven  by  an  appropriate  example  which 
Includes  the  selection  of  the  elgenstructure,  application  of 
elgenstructure  assignment,  and  evaluation  of  the  resulting 
closed-loop  system  by  pilots. 


Summary  of  current  Knowledge 

It  is  well  known  that  a  dynamic  system  and  in 
particular,  an  aircraft's  equations-of-motion,  can  be 
represented  by  n  first-order,  linear,  time-invariant, 
differential  equations  of  the  form: 


♦ 

X ( t )  =  Ax(t)  +  Bu( t )  (1.1) 

where 

x  =  nxl  state  vector 
u  =  mxl  control  vector 
A  =  nxn  plant  matrix 
B  =  nxm  input  matrix 

For  the  longitudinal  axis  of  an  aircraft,  the  A  and  B 
matrices  are  found  by  linearizing  the  vehicle's  equat lons-of- 
motlon  about  a  trim  condition  and  then  writing  the  stability 
axis  equations  in  matrix  form.  McRuer,  et  al.  contains  a 
rigorous  treatment  of  aircraft  equations-of-motion 
[5:Chap  5] . 

When  full-state  feedback  is  applied  to  a  system,  the 
knowledge  of  the  state  of  the  system  is  used  to  change  the 
input  in  order  to  make  the  output  behave  in  some  desired 
manner.  An  example  of  a  common  feedback  control  system  is 
the  cruise  control  available  on  cars.  In  this  application 
the  speed  of  the  car  is  the  output  of  the  system  and  is  fed 
back  to  the  input,  the  accelerator. 

In  multivariable  systems,  a  feedback  matrix  K  is 
formed  and  used  to  modify  the  control  vector,  u(t)  *  Kx(t) 

The  state  space  representation  can  now  be  written  as: 


The  quantity  (A  *  BK)  represents  the  closed-loop  system, 
and  the  problem  remains  to  determine  K  so  that  the  system 
has  desirable  behavior. 

As  is  well  known,  the  dynamic  response  of  a  linear 
system  is  determined  by  the  system's  eigenvalues  and 
eigenvectors.  Elbert  £61  or  Reid  [7]  contain  a  thorough 
treatment  of  linear  systems  theory.  The  eigenvalues  are 
simply  the  roots  of  the  characteristic  equation  for  the 
system: 

det ( A  +  BK  -  XI )  =  0  (1.3) 

i 

where  I  is  the  Identity  matrix  and  X  are  the  n 
eigenvalues.  The  eigenvalues  control  the  rate  of  decay  or 
growth  of  the  characteristic  modes  of  response  of  the  system. 
The  eigenvectors  are  found  by  solving: 

(A  +  BK)  p,  =  p.  X.  (1.4) 

where  pt  is  the  eigenvector  associated  with  the  1th 
eigenvalue.  The  eigenvectors  control  the  shape  of  the 
response  for  a  given  mode.  The  totality  of  the  eigenvalues 
and  eigenvectors  is  the  eigenstructure  of  a  system.  By 
picking  the  eigenstructure  for  a  system,  one  can  control  how 
the  system  will  respond  to  inputs,  disturbances,  initial 
conditions,  etc. 

Eigenvalue  placement  has  been  used  since  the  early 
1960's,  but  placing  the  eigenvectors  was  not  addressed  until 
1976  when  B.  C.  Moore  published  a  method  for  placing  both  the 
eigenvalues  and  eigenvectors  (81.  He  identified  the 
necessary  and  sufficient  conditions  for  simultaneous 
eigenvalue  and  eigenvector  assignment.  This  early  technique 
was  restrictive  in  that  it  did  not  accommodate  repeated 
eigenvalues,  and  it  permitted  only  certain  eigenvectors  to  be 


placed.  This  restriction  was  characterized  by  G.  Klein  as  a 
"class  of  generalized  eigenvector  chains"  [9:1401. 

Around  this  time,  many  control  system  engineers  realized 
that  while  full-state  feedback  permitted  a  good  deal  of 
flexibility  in  placing  the  elgenstructure,  the  necessity  of 
measuring  and  feeding  back  the  entire  state  was  not 
practical.  The  types  of  measurement  instrumentation  needed 
for  feedback  are  expensive,  sensitive  to  disturbances,  and 
require  frequent  calibration,  making  them  prohibitive  to  use 
in  the  flight  environment.  Since  the  outputs  are  already  a 
measured  variable  (all  systems  have  some  output  or  there 
would.be  no  point  to  the  system),  there  would  be  no 
additional  measurements  needed  for  the  feedback  signals. 
However,  since  the  number  of  outputs  is  usually  less  than  the 
number  of  states,  the  ability  to  control  the  elgenstructure 
placement  is  diminished. 

Although  this  study  used  full-state  feedback,  it  is 
interesting  to  briefly  examine  the  emergence  of  output 
feedback.  In  1978,  S.  Srinathkumar  wrote  the  benchmark  paper 
on  elgenstructure  assignment  using  output  feedback.  He 
defined  the  sufficient  conditions  under  which  the  eigenvalues 
and  eigenvectors  can  be  assigned  and  set  the  stage  for  the 
second  leap  in  the  use  of  this  design  technique  [10] . 
Srinathkumar  determined  the  maximum  number  of  eigenvalues  and 
eigenvectors  which  can  be  placed  for  a  given  system.  He 
concluded  that  "in  addition  to  assigning  min(n,  m+r-1) 
eigenvalues  (r-1)  eigenvectors  can  be  partially  assigned  with 
m  entries  in  each  vector  arbitrarily  chosen"  110:801.  In 
Sr inathkumar 's  paper,  m  is  the  number  of  inputs  and  r  is  the 
number  of  outputs.  For  the  state  space  model  shown  in 
Equation  (1.1),  Sr Inathkumar ' s  argument  would  correspond  to 
assigning  n  eigenvalues  and  a  maximum  of  m  elements  of  each 


eigenvector.  This  number  o£  eigenvalues  and  eigenvector 
elements  assumes  that  the  system  is  completely  controllable, 
which  is  not  always  the  case. 

The  Advanced  Fighter  Technology  Integration  (AFTI)  F-16 
aircraft  was  repeatedly  used  to  demonstrate  the  advantages  of 
output  feedback.  This  aircraft  has  additional  control 
surfaces  which  allow  it  to  fly  unlike  any  previous  airplane. 
The  additional  control  surfaces  mean  that  the  AFTI  system's 
control  matrix  is  of  higher  dimension  than  that  for  a 
conventional  aircraft.  For  this  design,  control  engineers 
had  multiple  inputs  Which  would  permit  enhanced  placement  of 
the  elgenstructure .  Designing  control  laws  for  the  AFTI  F-16 
became  very  popular  as  a  means  to  display  new  thoughts  and  • 
methods  on  output  feedback  for  elgenstructure  assignment. 

The  designs  successfully  decoupled  the  characteristic  modes 
of  the  aircraft  and  gave  the  engineers  insight  as  to  how  the 
eigenvectors  shaped  the  modes  of  a  system  [11,121. 

A  landmark  paper  was  written  by  Andry,  et  al.  in  1983 
which  summarized  elgenstructure  assignment  including 
full-state  and  output  feedback.  An  examination  was  also  made 
into  how  the  eigenvalues  and  eigenvectors  affected  the 
handling  qualities  of  the  aircraft;  that  is,  how  easy  or 
difficult  is  it  for  the  pilot  to  fly  the  aircraft  under 
various  conditions?  They  presented  eigenvectors  which  would 
generally  yield  good  handling  qualities  [13:7191.  Still, 
many  of  the  components  of  the  eigenvectors  were  assigned 
arbitrarily  and  much  of  the  eigenvector  placement  remains 
guesswork  or  based  on  experience. 

There  remains  a  good  deal  of  research  to  be  performed  in 
defining  the  relationship  between  elgenstructure  placement 
and  aircraft  handling  qualities.  Part  of  the  difficulty  in 
correlating  aircraft  handling  qualities  to  flight  control 


system  design  is  that  handling  quality  ratings  are,  for  the 
most  part,  the  subjective  evaluations  of  pilots.  Therefore, 
one  must  always  balance  a  pilot's  ratings  with  his  background 
and  flying  technique. 

The  specification  document  for  aircraft  handling 
qualities  is  written  in  parameters  which  are  easily 
interpreted  in  classical  control  theory  but  which  have  less 
direct  meaning  in  modern  control  theory.  Further  research 
into  this  relationship  and  more  specifically,  defining 
handling  qualities  parameters  in  terms  of  the  elgenstructure, 
would  make  it  much  easier  to  apply  this  technique  to  aircraft 
control  system  design.  Work  by  Stein  and  Henke  [14]  proposed 
a  method  for  forming  an  optimal  plant  matrix  which  could 
subsequently  be  decomposed  into  its  eigenvalues  and 
eigenvectors.  This  elgenstructure  could  then  be  used  in  the 
elgenstructure  assignment  procedure.  This  procedure  will  be 
examined  further  in  Section  III. 

Scope 

This  study  will  accomplish  each  key  step  in  the  flight 
control  design  procedure  in  order  to  prove  that  the  proposed 
elgenstructure  selection  process  is  applicable.  The  X-29A 
Advanced  Technology  Demonstrator  was  chosen  as  the  study 
aircraft  primarily  because  It  is  an  aerodynamically  unstable 
aircraft  that  requires  stability  augmentation,  and  also 
because  it  has  three  pitch  control  surfaces  which  will  allow 
fuller  assignment  of  eigenvector  elements. 

Extensive  use  will  be  made  of  the  Flight  Control 
Development  Laboratory's  simulation  of  the  X-29A.  This 
facility  Is  located  within  the  Flight  Dynamics  Laboratory  at 
Wr lght-Patterson  Air  Force  Base.  The  simulation  Is  a 
high-fidelity,  6  degree-of -freedom  simulation  which  drives 


the  Large  Amplitude  Multi-Mode  Aerospace  Research  Simulator 
(LAMARS)  motion-base  simulator.  This  simulation  was  used 
extensively  during  the  design  and  check-out  of  the  X-29A 
prior  to  its  first  flight  in  1984  and  was  used  more  recently 
to  examine  the  X-29A's  high  angle-of-attack  capabilities. 

The  simulation  incorporates  the  Grumman  aerodynamic  data  and 
flight  control  systems  and  includes  fourth-order  actuator 
models  on  an  analog  computer,  sensor  models  are  not  used  in 
this  simulation.  The  equations-of-motion  are  nonlinear,  and 
the  simulation  is  valid  in  the  entire  X-29A  flight  envelope. 

The  longitudinal  control  law  designed  in  this  study  will 
be  used  to  replace  the  Analog  Reversion  mode  longitudinal 
control  path.  The  Normal.  Digital  mode  was  not  selected  to  be 
replaced  due  to  the  increased  complexity  of  this  control  law 
which  Includes  fault  protection,  automatic  wing  camber 
control,  degraded  modes,  sensor  failure  contingencies,  and 
other  elements  which  must  be  taken  into  consideration  when 
developing  a  flight-ready  control  system. 

Full-state  feedback  will  be  used  in  the  control  system 
design.  This  method  may  not  be  practical  for  a  flight-ready 
application  because  of  the  expense  and  difficulty  of 
measuring  the  states,  however,  it  is  acceptable  here  since 
this  study  emphasizes  the  e lgenstructure  selection  process 
and  not  how  the  feedback  paths  would  be  implemented  on  an 
actual  aircraft. 

A  limited  number  of  up-and-away  flight  conditions  will 
be  examined  using  an  air-to-air  pitch  tracking  task  for  the 
piloted  portion  of  the  evaluation.  No  sensitivity  or  surface 
failure  analyses  will  be  done.  The  center-stick  feel  system 
characteristics  will  not  be  altered  from  those  of  the  actual 


The  development  and  evaluation  of  the  control  law  for 
this  study  can  be  separated  Into  five  distinct  steps: 

1.  Linearization  of  the  X-29A  airframe  and 
the  generation  of  the  state  space  matrices 

2.  Selection  of  the  desired  eigenstructure 
by  forming  an  ideal  plant  matrix  and 
performing  a  spectral  decomposition  of 
this  matrix 

3.  utilization  of  eigenstructure  assignment 
to  obtain  the  feedback  gain  matrix 

4.  Generation  of  a  pilot  model  to  accomplish 
an  off-line  analysis  of  the  closed-loop 
system 

5.  '  Implementation  of  the  control  law  on 

the  LAMARS  and  piloted  evaluations  of 
the  aircraft 

Each  of  these  steps  will  be  discussed  briefly  in  the 
remainder  of  this  section. 

Linearization.  A  linearized  model  of  the  X-29A  airframe 
was  obtained.  In  order  to  Include  both  the  rigid  body 
derivatives  and  flexible  body  increments  In  the  linear  model, 
and  to  have  an  accurate  representation  of  the  LAMARS 
simulation,  the  simulations 's  aerodynamic  model  was 
linearized  rather  than  using  existing  Grumman  or  NASA 
linearized  aerodynamic  data.  The  state  variables  and  control 
surfaces  were  perturbated  about  trim  equilibrium  conditions, 
and  nondlmenslonal  stability  and  control  derivatives  were 
obtained.  The  longitudinal  state  variables  used  for  this 
study  were : 

u  -  forward  velocity  (feet  per  second) 

a  -  angle-of-attack  (radians) 

q  -  pitch  rate  (radians  per  second) 

9  -  pitch  angle  (radians) 


tamtam 


The  three  pitch  control  surfaces  are  the  close-coupled, 
all-moving  canard;  variable  camber,  tral 1 lng-edge  flaps;  and 
the  strake  flaps.  The  nondimensional  stability  and  control 
derivatives  were  transformed  into  the  dimensional  derivatives 
used  to  calculate  the  elements  of  the  A  and  B  matrices  of 
the  linear,  time-invariant  system. 

Elgenstructure  Selection .  Choosing  the  desired 
eigenvalues  and  eigenvectors  of  the  augmented  matrix  was  one 
of  the  most  challenging  parts  of  this  study.  Sufficient 
empirical  data  exists  to  choose  eigenvalues  which  should 
yield  good  response  decay  rates.  However,  when  picking  the 
eigenvectors,  in  the  past,  one  relied  upon  a  few  broad 
statements  (such  as  orthogonality  of  the  short  period  and 
phugoid  modes,  and  the  null  participation  of  certain  vectors 
in  certain  responses)  and  experience.  This  study  generated 
the  eigenvalues  and  eigenvectors  by  building  an  optimal  A 
matrix  based  on  knowledge  about  the  relationship  between 
aircraft  response  performance  and  the  dimensional  stability 
derivatives.  This  ideal  matrix  was  then  decomposed  into  its 
eigenvalues  and  eigenvectors,  and  these  became  the  desired 
elgenstructure  of  the  control  law  design. 

Control  Law  Design .  Elgenstructure  assignment  using 
full-state  feedback  was  used  to  develop  a  feedback  gain 
matrix,  K  .  A  computer  program  called  MODES  (see  Ref.  15 
for  a  description)  was  used  to  perform  this  task  in  a  quick 
and  efficient  manner.  Since  all  the  states  were  available  in 
the  simulation,  full-state  feedback  was  easily  accomplished. 
The  closed-loop  system  matrix,  (A  +  BK)  ,  decomposed  into  the 
desired  elgenstructure  described  in  the  previous  section.  In 
this  case,  the  A  and  B  matrices  were  the  bare  airframe 
X-29A  plant  and  control  matrices,  respectively. 


Off-line  Analysis .  A  linear  pilot  model  was  used  to 
check  that  the  pilot-aircraft  dynamic  system  exhibited  good 
performance.  The  pilot  model  used  was  sufficient  to  account 
for  the  capabilities  and  limitations  of  the  pilot  In  the 
closed-loop  system.  The  form  of  the  pilot  model  used  In  this 
study  Is  based  on  the  Neal-Smlth  report  [161: 

K  e_xs(T  s  +  1) 

p  p . 

Y  =  - 1 - 

P  {  T  S  +  1  ) 

P  2 

where 

Yp  =  pilot  model  transfer  function 

Kp  =  pilot  gain 

e~Ts  =  pilot's  neuromuscular  time  delay,  t 

t  =  pilot  lead  added 

pi 

=  pilot  lag  added 

2 

LAMARS  Simulation.  The  final  version  of  the  control  law 
was  Integrated  into  the  LAMARS  simulation.  Once  the  Analog 
Reversion  control  system  had  been  modified  with  the  new 
longitudinal  control  law,  digitally  controlled  Inputs  into 
the  control  system  were  used  to  generate  time  histories. 

These  were  used  to  check  for  proper  implementation  of  the 
control  law  and  for  any  additional  time  delay  which  may  have 
been  introduced.  An  air-to-air  tracking  task  was  employed  to 
obtain  Cooper-Harper  handling  quality  ratings  from  the  pilot 
during  the  man-ln-the-loop  simulation  runs.  The 
Cooper-Harper  scale  Is  an  accepted  pilot  ratings  tool  that 
has  been  used  In  flight  testing  for  over  15  years  [17].  A 
correlation  between  Cooper-Harper  ratings  and  flying  quality 
level  is  addressed  in  MIL-F-8785 C.  A  sufficiently  large  test 
matrix  was  used  in  the  piloted  sumulations  in  order  to  obtain 


(1.5) 


=  0.3  sec 


The  primary  resources  necessary  to  accomplish  this  study 
were  contained  within  the  LAMARS  simulation  facility.  The 
X-29A  simulation  uses  a  Gould/SEL  32/77  and  32/97  digital 
computer,  as  well  as  an  EAI  781  analog  computer  connected  by 
an  EAI  Hyshare  system.  A  CSPI  MAP  multi-array  processor  is 
used  for  the  aerodynamic  table  look-ups.  The  LAMARS  itself 
is  a  five  degree-of-f reedom,  beam-type  motion  simulator.  The 
20-foot-dlameter  dome  contains  a  single-seat  fighter  cockpit 
and  spherical  dome  display  with  a  target  projector.  The 
cockpit  was  configured'  as  a  generic  fighter  for  this  study 
with  a  head-up  display,  center  stick,  throttle  quadrant,  and 
standard  front-panel  Instrumentation.  The  center-stick  was  * 
programmed  to  have  the  X-29A's  feel  system  characteristics 
such  as  the  stick  force  gradient  and  break-out  forces.  The 
head-up  display  presented  a  fixed  reticle  for  the  pilot  as 
well  as  the  pitch  command  tracking  bar. 

The  interactive  design  program  TOTAL,  hosted  on  the  ASD 
Cyber  system,  was  used  to  manipulate  transfer  functions  and 
obtain  preliminary  Nlchol's  plots.  Also,  as  previously 
mentioned,  the  MODES  program  was  used  to  find  the  feedback 
matrix.  This  program  is  available  either  on  an  AFIT 
minicomputer,  or  it  can  be  used  on  a  personal  computer. 


Linearization  of  the  X-29A 


II . 

One  of  the  basic  elements  needed  to  perform 
elgenstructure  assignment  is  a  mathematical  model  of  the 
unaugmented  system  in  state-space  form.  Since  pilot  handling 
quality  ratings  were  to  be  taken  from  the  LAMARS 
man-in-the-loop  simulation  of  the  X-29A  supplemented  with  the 
control  law  designed  in  the  course  of  this  project,  it  was 
necessary  to  obtain  an  accurate  model  of  the  LAMARS 
simulation,  rather  than  using  existing  NASA  or  Grumman 
linearized  data.  This  was  necessary  since  models  of  the 
X-29A  aerodynamics  used  by  the  various  agencies  varied  in  the 
level  of  sophistication  and  the  revision  level  implemented. 
The  aerodynamic  model  used  on  the  LAMARS  is  referred  to  as 
AER07B  by  the  X-29A  community.  It  is  a  six  degree-of-f reedom 
nonlinear  model  which  Includes  flexibility  Increments  and 
contains  over  60,000  data  points.  It  is  valid  in  the  entire 
X-29A  envelope  and  is  modified  to  handle  high  aircraft 
angles-of-attack  (-50  to  +  90  degrees  AO A) .  For  this  study, 
four  flight  conditions  were  examined  and  feedback  gain 
matrices  were  found  for  each.  The  elements  of  the  gain 
matrices  were  then  curve  fit  between  the  flight  conditions  to 
provide  a  piecewise  linearization.  The  flight  conditions 
studied  were  Mach  =  0.6,  0.8  and  altitude  =  10K,  30K  feet. 

Formulation  of  the  Linearized  Equations  of  Motion 

The  linearized,  perturbation  equat ions-of-motion  were 
taken  from  Reference  5.  The  equatlons-of-motion  were  needed 
in  the  body  reference  frame  for  the  X-29  linearization  but 
also  were  needed  in  the  stability  axis  frame  for  use  in  the 
development  of  the  optimal  elgenstructure,  discussed  in 
Section  III.  Only  the  body  axis  equations  will  be  discussed 


here  and  the  body  axis  frame  Is  shown  In  Figure  1.  A  number 
of  assumptions  were  made  during  the  formulation  of  the 
equations  and  are  listed  below: 

1)  The  airplane  Is  a  rigid  body. 

2)  The  earth  is  fixed  in  apace. 

3)  The  airplane  has  constant  mass  and  mass  distribution. 

4)  The  XZ  plane  is  a  plane  of  symmetry. 

5)  All  angles  are  small. 

6)  Steady-state  lateral  trim  conditions  are  P,  =  R,  =  V, 

=  *,  =  0. 

7)  The  flow  Is  quasisteady. 

8)  Variations  of  atmospheric  properties  are  negligible 
for  small  altitude  perturbations. 

9)  The  earth  Is  flat. 

The  resulting  equations,  in  the  body  axis  frame,  are  shown 
below: 

u  -  a  =  Xu  u  +  Xa  ct  +  (Xq  -  WJq  -  g  cos  8 ,  8 

+  Xfig6s  +  X^5r  +  Xs_fi.  (2.1 

<u.  -  Z£)a  =  Zu  u  +  Za  a  +  (Zq  +  U. )q  -  g  sin  8,  6 

+  ZS  £S  +  h  Sr  +  ZS  Sc  (2-2 

s  r  c 

-M*  i  +  q  =  M  u  +  M..  ot  +  M_  q  +  M-  o  „  +  Mr  c  „  +  M*  o ,  (2.3 

—  u  u-  <3  °ss  ^rr  c 

6  =  q  (2.4 

In  these  equations,  6,5,  and  6.  are  the  strake  flap, 
symmetric  flaperon,  and  canard  surface  deflections, 
respectively.  These  control  surfaces  are  shown  In  Figure  2. 

The  equations-of-motion  were  needed  in  state-space  form 
for  the  e lgenstructure  assignment  process.  This  was  quickly 
accomplished  and  Equation  (2.5)  is  the  resulting  form. 


-29A  Control  Surface  Definition 


where 


X 


A  computer  program  was  written  which  quickly  and 
efficiently  computed  the  dimensional  derivatives  and  formed 
the  A  and  B  matrices  for  the  LAMARS  X-29A  simulation. 


First  the  aircraft  was  trimmed  in  one  of  the  four  flight 

conditions  of  Interest.  In  the  analog  reversion  mode,  the 

strakes  and  canards  are  set  to  zero  and  only  the  flaperons 

are  used  to  balance  the  forces  when  trimming.  From  the 

trimmed  condition,  one  of  the  states  was  changed  first 

through  a  small  positive  perturbation,  then  through  an  equal 

negative  perturbation  and  the  total  aerodynamic 

nondlmenslonal  coefficients  (Cx,  Cz,  CM)  were  determined  from 

the  X-29A  simulation's  aerodynamic  model  at  each  perturbation 

step.  Once  this  was  accomplished,  the  nondlmenslonal 

stability  derivatives  were  calculated  (Cv  ,  C_  ,  C_  ,  etc) 

xct  zu  mq 

based  on  the  difference  between  the  aerodynamic  coefficients 
for  the  positive  and  negative  perturbations  and  the  size  of 
the  perturbation.  The  variables  perturbated  were:  alpha, 
alpha-dot,  pitch-rate,  Mach  number,  strake,  flaperon,  and 
canard.  It  should  be  noted  that  these  coefficients  were 
measured  In  the  body  axis  frame. 

Next  the  dimensional  derivatives  were  formed  from  the 
nondlmenslonal  derivatives  using  the  following  equations: 


psu.  M  w„ 

'  7  V  C*  *  V 

1/sec 

(2.6) 

psu:  w.  w 

xa  =  TnT  «V  2*r(CX  +  ?  CxMn 

ft/sec2 

(2.7) 
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where  VT.  is  the  total  velocity  in  ft/sec 

And  similarly  for  the  Z  and  M  derivatives. 

With  the  dimensional  derivatives  in  hand,  the  A  and  B 
matrices  were  formed  using  Equation  (2.5).  This  process  was 
repeated  for  each  flight  condition.  The  A  and  B  matrices 
for  Mach  =  0.6  and  altitude  =  10k  feet  are  shown  below: 


A  = 


B  = 


- .  1262E-01 

.  3680E+02 

-. 3451E+02 

-. 3213E+02 

- .  8438E-04 

- . 1048E+01 

. 9906E+00 

- . 2657E-02 

- .  5757E-03 

. 1687E+02 

-. 4844E+00 

. 1553E-03 

.  0000E+00 

. OOOOE+OO 

. 1000E+01 

.  0000E+00 

'  . 1590E+01 

. 5594E+01 

- . 7027E+01 

-.  4534E-01 

- .  3359E+00 

- .  6440E-01 

- .  2990E  +  01 

- . 5455E+01 

.  8449E+01 

. 0000E+00 

.0000E+00 

.0000E+00 

The  computer  program  also  accessed  a  subroutine  which 
determined  the  eigenvalues  of  the  A  matrix.  The  dynamic 
characteristics  for  the  character istic  modes  were  calculated 
from  the  eigenvalues.  For  the  above  flight  condition  ,  the 
eigenvalues  and  dynamic  characteristics  for  the  bare  airframe 
are 


h  = 


3.334  +  j(0.0) 

X2  =  -4.865  +  j(0.0) 

X3  4  =  -.007138  ±  J(0. 06258) 

C  =  .11332 


Unstable  Short  Period 
Phugoid 


w  =  .06299  rad/sec 


While  it  has  taken  only  a  short  section  to  describe  the 
linearization  process  used  to  obtain  the  unaugmented  A  and 
B  matrices  of  the  X-29A,  it  should  be  noted  that  this  was  a 
very  critical  and  laborious  portion  of  the  project.  The 
importance  of  having  an  accurate  state-space  model  has 
already  been  discussed. 

The  accuracy  of  the  linearization  process  was  checked 
by  comparing  the  elements  of  the  A  and  B  matrices 
developed  here  to  those  appearing  in  a  NASA  X-29A  linear 
analysis  report  [18].  The  correlation  was  very  good  and  the 
next  step  In  the  project  was  to  develop  a  method  to  select 
desired  eigenvalues  and  eigenvectors  for  the  closed-loop 
system. 
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III.  Development  of  optimal  Elgenstructure 
Background 

As  has  been  previously  argued  In  Section  I ,  there  is 
little  In  the  current  handling  quality  specification  which 
can  guide  the  control  system  designer  wanting  to  use 
elgenstructure  assignment.  Desirable  eigenvalues  can  be 
computed  based  on  specifications  of  the  system  damping  and 
natural  frequency,  however  there  is  no  Information  In  the 
specifications  which  can  be  used  to  directly  compute  the 
desired  eigenvectors.  One  of  the  major  challenges  of  this 
effort  was  to  derive  a  desirable  elgenstructure  from  the 
handling  quality  specif ications . 

Stein  and  Henke  proposed  a  design  procedure  using 
quadratic  optimal  control  theory  and  a  handling  quality 
oriented  cost  function  for  flight  control  systems.  They  used 
data  from  MIL-F-8785  to  define  a  system  of  differential 
equations  having  the  same  form  and  desirable  handling 
qualities  of  the  augmented  aircraft  [14:  Sec  III.  The 
portion  of  their  work  of  Interest  here  Is  how  they  formed  a 
plant  matrix  possessing  desirable  handling  qualities.  Stein 
and  Henke  applied  their  design  procedure  to  the  lateral- 
directional  axis.  The  following  discussion  extends  their 
work  to  the  longitudinal  axis,  with  the  intent  of  using 
full-state  feedback. 

Handling  Quality  Model 

The  beginning  point  for  generating  the  handling  quality 
model  is  the  A  (plant)  matrix  for  the  rigid-body,  longi¬ 
tudinal  equatlons-of-motion,  written  in  the  stability  axis 
system.  We  want  an  A  matrix  which,  when  placed  in  the 


x  =  A  x  system,  will  yield  desired  handling  qualities  as 
defined  by  the  current  MIL-PRIME  Standard.  This  A  matrix 
will  provide  the  desired  e igenstructure  which  will  be  used  in 
the  full-state  e igenstructure  assignment. 

Since  the  elements  of  the  A  matrix  are  composed  of  the 
dimensional  derivatives,  the  problem  of  forming  the  optimal 
matrix  reduces  to  assigning  values  to  these  derivatives  such 
that  the  optimal  matrix  results.  If  the  dimensional 
derivatives  can  be  related  to  handling  quality 
specifications,  then  the  optimal  values  of  the  dimensional 
derivatives  can  be  solved  for. 

There  are  twelve  longitudinal,  aerodynamic  dimensional 
derivatives.  In  order  to  solve  explicitly  for  these  twelve 
variables,  one  would  need  twelve  linearly  independent 
equations.  Fortunately,  four  of  the  derivatives  can  be 
assummed  zero : 

-  is  primarily  due  to  the  change  in  drag  on  the  pitch 
control  surface  which  is  small  in  comparison  to  the 
total  drag  15:273] 

-  is  due  to  the  change  in  lift  on  the  pitch  control 
surface  primarily  as  a  result  of  an  aerodynamic 
time  lag  [5:274] 

Xq  -  again,  the  increased  drag  is  negligible  to  the 
first  approximation  [5:275] 

Zq  -  again,  the  increased  lift  has  little  effect  and 
1  is  neglected  [5:277] 

This  leaves  eight  dimensional  derivatives  which  must  be 
solved  for  from  eight  as  yet  unknown  equations.  The  approach 
taken  to  identify  the  eight  equations  was  to  sift  through  the 
MIL-PRIME  Handbook  and  other  handling  quality  sources  and 
identify  specifications  which  could  be  related  to  an 


m 


expression  involving  the  dimensional  derivatives.  After 
several  passes  through  the  documents,  eight  equations  were 
found  and  are  discussed  below,  along  with  the  handling 
quality  specification  which  corresponds  to  each  equation. 

When  examining  specif lcations  from  the  literature,  the  flight 
phase  considered  was  Category  A  (alr-to-alr  combat)  and  the 
aircraft  type  was  Class  IV  (hlgh-maneuverablllty). 

The  MIL-PRIME  Handbook  provides  recommended  values  for 
the  phugoid  mode  damping  ratio,  Cph  14:186).  McRuer  gives  an 
approximate  form  for  the  phugoid  damping  (5:336): 


s  -X  - 
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-  9) 

O.M,  - 


(3.1) 


The  MIL-PRIME  Handbook  doesn’t  address  the  natural 
frequency  of  the  phugoid,  but  a  review  of  various  aircraft's 
dynamic  characteristics  can  provide  good  values.  The 
approximate  form  for  the  phugoid  natural  frequency  Is  (5:336) 
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Recommended  values  for  the  short  period  damping  are 
found  In  MIL-F-8785C  (2:13)  and  the  short-period 
approximation  gives  (5:336) 

2^sp  s  +  +  <3’3> 

The  short-period  natural  frequency  is  partially 
restricted  by  allowable  values  for  the  Control  Anticipation 
Parameter  (CAP)  which  is  directly  proportional  to 
(4:1921.  Also  Etkln  gives  a  range  of  values  for  u,p  based  on 
experiments  using  a  variable  stability  aircraft  (19:513). 

The  short-period  approximation  for  is  (  5:336) 
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Allowable  values  for  the  normal  acceleration  change  per 
unit  angle-of -attack,  n/a,  are  given  In  MIL-F-8785C  [2:141. 
An  approxlnvate  form  for  n/ot  Is  given  by  Roskam  [  20:  5381  : 


n/a  ^  -2a/g 


(3.5 


The  modal  response  ratios  have  a  direct  bearing  on  the 
eigenvectors.  For  the  phugoid  mode,  Etkin  [19:3241  and 


McRuer  [5:3501  both  state  that  there  is  a  90*  phase 
difference  between  speed,  u,  and  the  pitch  angle,  0.  This 
angle  can  be  approximated  by  [5:3501 


Oox  +  (g  m  /m  , ) 


For  the  short-period  mode,  Etkin  [19:324]  shows  that  the 
phase  difference  between  a  and  8  Is  small  and  the  angle  can 
be  approximated  by  [5:347] 
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The  static  to  short-period  gain  ratio  must  be  greater 
than  allowable  limits,  as  Indicated  in  Reference  1,  page  66 
This  corresponds  to  defining  a  minimum  low-frequency  gain 
margin.  The  static  to  short-period  gain  ratio  can  be  be 
expressed  as  [5:4241 
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These  are  the  eight  equations  in  eight  unknowns  that  are 


required.  When 


or  “oh  appeared  on  the 


right-hand  aide  of  an  equation,  they  were  replaced  by  the 
appropriate  approximation  expression. 

Three  sets  of  data  were  chosen  for  the  handling  quality 
parameters  (the  left-hand  side  of  the  equations):  one  set  to 
yield  good  handling  qualities,  one  to  yield  fair,  and  one  to 
yield  poor.  These  three  sets  correspond  to  Level  1,  2,  and  3 
respectively,  as  defined  in  the  MIL-PRIME  Handbook  (4:96). 
The  values  chosen  for  the  handling  quality  parameters  are 
shown  in  Table  I. 


Table  I 

Handling  Quality  Parameter  Values 


Good 

Fair 

Poor 

mam 

0.007 

0.002 

-0.004 

U)2  , 
p  fc 

0.0025 

0.0025 

0.0025 

2tw 

_ '  1  p _ 

4.90 

1.50 

0.40 

■EM 

12.25 

4.0 

1.0 

n/ct 

30.00 

20.00 

20.0 

B 

85.00 

70.00 

45.0 

QBS 

15.00 

24.00 

45.0 

3-0 

— 

1.0 

0.01 

2C  w  .  &  ui ‘  v 

*  p  n  p  n 


A  value  of  0.05  was  chosen  for  the 


phugold  natural  frequency  based  on  a  survey  of  fighter 
aircraft.  This  value  was  used  in  all  three  sets  of  data  in 
order  to  Isolate  the  effect  of  changing  the  damping  ratio  in 
the  phugold  mode.  The  MIL-PRIME  directive  was  consulted  in 
picking  values  for  of  0.07,  0.02,  and  -0.04.  Note  that 
the  last  value  will  produce  an  unstable  phugold  with  a 
tlme-to-double  of  347  secor/is,  which  falls  within  the  Level  3 
criteria.  Since  the  phugold  mode  characteristics  are 


unlikely  to  be  noticed  by  the  pilot  In  an  air-to-air  tracking 
task,  the  three  different  sets  of  values  may  have  no  effect 
on  pilot  comments  and  handling  quality  ratings. 

2Cuis  &  ■,/  p  .  MIL-F-8785C  gives  a  range  of  acceptable 

values  for  the  short-period  damping  ratio.  The  Level  1  value 
picked  was  0.7,  Level  2  was  0.375  (this  also  falls  within  the 
lower  range  of  acceptable  Level  1  values  but  was  felt  to  be 
low  enough  that  when  coupled  with  other  Level  2  parameters 
would  not  give  Level  1  handling  qualities),  and  Level  3  was 
0.2.  The  values  for  the  short-period  undamped  natural 
frequency  were  selected  to  yield  appropriate  values  for 
Control  Anticipation  Parameter  (CAP),  which  range  between 
0.28  and  1.0  rad/g*s;  for  fine  tracking.  The  values  for  u)sp 

were  3.5  for  Level  1,  2.0  for  Level  2,  and  1.0  for  Level  3. 

n /a  .  Values  were  selected  from  Figure  13  of  the 

MIL-PRIME  Handbook  using  the  u  selected  for  each  level. 

s  p 

Modal  Phase  Angles  .  The  angles  were  chosen  based  on 
generalizations  about  how  the  vectors  in  the  Argand  diagram 
behave  in  the  short-period  and  phugoid  modes.  Since  these 
angles  have  a  direct  influence  on  the  eigenvectors,  the 
choices  for  the  angles  are  likely  to  have  a  strong  effect 
during  the  eigenstructure  assignment. 

Static  to  Short-Period  Gain  Ratio  .  Values  were 
selected  to  give  9  db  of  low-frequency  gain  for  Level  1,  0  db 
for  Level  2,  and  -40  db  for  Level  3. 

Simultaneous  Solution  for  HQ  Mode  1 

The  eight  non-linear  equations  were  solved  simultaneous¬ 
ly  using  a  routine  from  the  International  Mathematical  « 
Statistics  Libraries  (IMSL)  software  package.  The  routine 
uses  the  secant  method  for  finding  the  roots  to  an  equation 
and  requires  the  user  to  provide  an  initial  guess  for  each  of 


the  design  variables  (the  dimensional  derivatives  for  this 
application).  Choosing  the  initial  conditions  proved  to  be  a 
difficult  task  when  trying  to  solve  all  eight  equations 
simultaneously.  If  the  algorithm  converged  at  all,  the  roots 
were  usually  unreasonably  large  (the  algorithm  converged  on 
an  unacceptable  local  minimum).  Often  the  algorithm  Iterated 
a  design  variable  to  the  wrong  sign  (for  Instance  a  positive 
Ma)  and  an  error  would  result  for  trying  to  take  the  square 
root  of  a  negative  number.  Due  to  these  difficulties,  it  was 
decided  to  reduce  the  number  of  equations  being  solved 
simultaneously  and,  as  valid  solutions  were  found,  to  add 
back  equations  one  or  two  at  a  time  until  all  the  equations 
were  being  solved  together. 

Using  the  solution  for  a  given  dimensional  derivative  as 
the  initial  guess  for  that  design  variable  the  next 
iteration,  allowed  swift  and  valid  convergence  of  the 
algorithm.  Unfortunately,  there  is  not  a  subset  of  the  eight 
equations  which  would  give  two  equations  in  two  unknowns  (or 
three  in  three).  Therefore,  Equations  (3.3),  (3.4),  and 
(3.7)  were  chosen  to  be  solved  simultaneously  and  values  for 
M*  were  supplied  in  order  that  only  three  unknowns  remained 
(H,,  M  ,  &  Z „ . )  .  M‘  was  chosen  to  be  given  a  constant  value 
because,  as  indicated  in  Table  5-6  of  Reference  5,  it  has  the 
least  widespread  effect  on  the  dynamics.  The  three  equations 
were  solved  several  times  with  different  values  of  to 
insure  that  the  solutions  found  weren't  overly  sensitive  to 
variations  in  M.‘, . 

J.- 

The  values  found  for  Ma,  ,  and  Za  were  used  as  the 
initial  guess  when  another  equation  (3.6)  was  added  to  the 
system  (adding  only  one  more  unknown).  This  process  was 
continued  until  values  were  found  for  all  eight  variables  at 
four  separate  flight  conditions  and  for  each  HQ  level.  In 


the  end,  only  seven  of  the  equations  were  solved  simultane¬ 
ously;  Equation  (3.5)  was  not  used.  M^  was  not  allowed  to 
optimize  for  the  reasons  stated  earlier.  The  four  flight 
conditions  examined  were  Mach  =  0.6,  0.8  and  altitude  =  10K, 
30K.  The  results  for  the  dimensional  derivatives  are  shown 
In  Table  II.  It  Is  pleasing  to  see  that  some  of  the 
dimensional  derivatives  are  nearly  constant  with  flight 
conditions.  Hopefully,  this  will  lead  to  a  feedback  gain 
matrix  which  will  require  limited  gain  scheduling  with  flight 
condition  (dynamic  pressure  typically). 


Optimal  A  Matrix  and  Elgenstructure 


With  the  necessary  dimensional  derivatives  in  hand  now, 
the  optimal  A  matrices  were  formed  for  each  of  the  twelve 
sets  of  data.  The  matrix  for  M  =  0.6,  Alt  =  10K  for  each 
handling  qualtity  level  is  shown  below  (in  stability  axis): 


Level  1 


f-  .  7200E-0  2  0 . 2200E+02  O.OOOOE+OO  -.3217E+02 

j -  .  1276E-03  - . 1315E+01  0.1000E+01  0.0000E+00 

!  -  .  2053E-05  - . 7537E+01  -.  3587E+01  0.0000E  +  00 

1 0 . 0000E  +  00  0.0000E  +  00  0.1000E  +  01  O.OOOOE+OO 


Level  2 


-  .  1000E-01  0 . 1605E  +  02  0.0000E  +  00  -.3217E+02 

- .  5972E-03  - . 9616E  +  00  0.1000E+01  0.0000E+00 

- .  1841E-02  - . 3482E  +  01  -.  5393E+00  0.0000E+00 

O.OOOOE+OO  O.OOOOE+OO  0.1000E+01  O.OOOOE+OO 


-J-  Aj.  r..  -J,  -J,  rj.'j'j.  v»V.  rAv 


Level  3 


- . 2130E-01  - . 1164E+02  O.OOOOE+OO  -.3217E+021 

- . 1868E-02  - . 1022E+01  0.1000E+01  O.OOOOE+Oo' 

-.2914E-02  -.163  5E  +  01  0.6208E+00  0.0000E  +  00 

0.0000E+00  O.OOOOE+OO  O.IOOOE+Ol  O.OOOOE+OO 

It  was  important  to  establish  that  these  matrices  did  indeed 
have  the  dynamic  characteristics  to  place  them  in  the 
handling  quality  level  for  which  they  were  designed.  For 
this  reason,  the  eigenvalues  were  found  for  each  matrix.  It 
was  verified  that  two  characteristic  modes  were  present  and 
the  damping  ratio  and  natural ; frequancy  were  calculated.  In 
addition,  several  other  handling  quality  criteria  were  used 
to  check  the  optimality  of  each  matrix: 

1)  Check  CAP  vs  ssp  from  Figure  3  [4:193] 

where:  CAP  =  (g  ♦  s  1  ) 

n/a  =  -Za/g  ( g«  rad' 1  ) 

2)  Check  wspT0  vs  Csp  from  Figure  4  [4:211,386] 

where:  T0  =  <3ec> 

3)  Check  w  vs  n/a  from  Figure  5  [4:214] 

*  p 

4)  Check  that  1/TP  is  greater  than  0.38  for  Level  1  and 

2 

0.24  for  Level  2  [4:385] 

5)  Check  that  g-jj  is  less  than  0.06  for  Level  1,  0.15  for 
Level  2,  and  0.24  for  Level  3 

where:  gj  =  |  [xu  -  (Xa  -  g)-^-]  (deg'knot'1  ) 

6)  Check  vs  C 

-  „  ,  •  5  P 


from  Figure  6  [18:514] 


The  dynamic  characteristics  o£  all  12  matrices  are  summarized 
In  Table  III . 


Most  of  the  data  falls  fairly  well  into  the  handling 
quality  level  for  which  its  matrix  was  designed.  The  ujpTQ 

vs  ?  check  gives  borderline  Level  1/Level  2  handling 
qualities  for  data  sets  5-8,  and  borderline  Level  2/Level  3 
for  data  sets  9-12.  Also,  for  all  cases  falls  within 
Level  1  criteria.  Again,  it  is  felt  that  this  criteria  will 
not  be  a  factor  in  the  air-to-air  tracking  task. 

The  damping  ratio  for  the  Level  3  data  sets  were 

i 

unanticipated.  By  the  phugoid  approximation  we  should  have  a 
damping  ratio  of  -0.04;  what  was  produced  was 
s  h  =  -0.9404.  Clearly  Equation  (3.1)  was  not  a  good 
approximation  in  this  case.  This  is  substantiated  by  the 
fact  that  the  condition  for  validity  for  the  approxlmat ion 
was  erroroneously  violated  [5:3361.  The  resulting 
t lme-to-double  of  15  seconds  is  too  fast  for  Level  3  and 


would  probably  be  considered  uncontrollable.  However,  it  was 
decided  that  it  would  be  interesting  to  leave  the  damping 
ratio  as  is  for  these  cases  to  see  if  the  pilot  would  notice 
the  poor  phugoid  damping. 

The  A  matrices  were  considered  to  be  optimal  at  this 
point  and  a  spectral  decomposition  was  performed  to  get  the 
eigenvalues  and  eigenvectors.  Again  for  Mach  =  0.6, 

Alt  =  10K,  the  e igenstructure  for  each  level  is  shown  In 
Table  IV.  The  next  step  in  the  control  system  design  process 
was  to  find  the  feedback  gain  matrix  given  the  bare  airframe 
X-29A  A  &  B  matrices,  and  the  just  obtained  optimal 
e igenstructure . 
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Table  III 

Optimal  Matrices'  Dynamic  Characteristics 


Table  IV 

Optimal  Eigenstructure 
(Mach  =  0.6,  Alt  =  10k 


Level 

Short-Period 

Phugoid  | 

= 

-2.45 

± 

2.4994 

X  h 

— 

-.0035 

± 

.0501 

4  t  p  : 

= 

-.925 

± 

-1.805 

*phS 

= 

-618.7 

307.8 

1 

;'ip: 

= 

0.084 

± 

-.8021 

^ph: 

= 

0.0232 

± 

-.012 

*,p3 

= 

1.909 

± 

1.1221 

<ph3 

= 

-.0480 

± 

0.251 

^  s  p  4 

= 

-.153 

± 

-.6138 

r 

^  p  h  4 

= 

0.5659 

± 

0.919 

% 

— 

-0.75 

± 

1.8573 

X  . 

- 

-.0015 

± 

.0498 

_£>  - 

•V 

s  s  p  i 

= 

-5.03 

± 

2.7050 

*pM 

= 

-852.8 

± 

-111. 

2 

■  spi 

s 

-0.48 

± 

0.4326 

*pht 

= 

0.4617 

.0577 

r 

S  _  _ 

= 

-0.91 

± 

-.7988 

^  p  h  } 

s 

-.0689 

± 

.0122 

<S04 

— 

-0.20 

0 . 5687 

^  0  h4 

= 

0.2853 

± 

1.377 

X 

_*  2. 

= 

-.257 

± 

0.9828 

X  k 

Eh  _ 

= 

0.0456 

± 

.0165 

*■ 

L  ;  p  1 

= 

-29.3 

i 

15.257 

^  p  h  1 

= 

-3996. 

± 

-807. 

3 

t' 

•  s  p  i- 

= 

-.323 

± 

0.5842 

= 

7.1925 

± 

1.509 

^  S  p ; 

= 

-.876 

i 

0.1583 

^Ph: 

0.1891 

± 

0.232 

r 

'  S  p  4 

= 

0.369 

± 

0.7948 

2“’ 

W'  p  h  4 

5.2952 

± 

3.170 

r»nWi 
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Eiqenstructure  Asslqnment 


Once  the  unaugmented  aircraft  model  and  the  optimal 
e igenstructure  were  found  in  Sections  II  and  III, 
respectively,  the  feedback  gain  matrices  could  be  calculated. 
This  was  accomplished  using  the  MODES  software  package  [151. 

Theory  Behind  MODES 

MODES  calculates  the  feedback  gain  matrix  K  which  will 
cause  the  closed-loop,  linear,  time-invariant  system 


x ( t )  =  l A( t )  +  B(t)K(t) ]x(t) 


(4.1) 


to  have  the  desired  elgenstructUre  selected  in  Section  III. 
Since  the  number  of  inputs  (3)  did  not  equal  the  number  of 
states  (4)  for  this  problem,  there  was  no  exact  solution  to 
the  problem 


(4.2) 


K  =  B''  [PoAPc1  -  AJ 

where 

Pj.  is  the  matrix  of  desired  eigenvectors 
A  is  the  diagonal  matrix  of  desired  eigenvalues 

and  A  =  A  +  BK  =  PDAP~‘ 


Therefore,  the  eigenvectors  picked  In  Section  III  were 
not  fully  realizable  and  some  projection  of  these 
eigenvectors  was  needed.  Singular  value  decomposition  was 
used  to  find  the  range  space  of  the  achievable  eigenvectors 
of  the  augmented  matrix  [A  -  X^l|Bj  ,  where  V  is  a  desired 
eigenvalue.  Once  the  range  space  was  found,  singular  value 
decomposition  was  used  again  to  find  the  projection  of  the 
desired  eigenvector  into  the  range  space  of  achievable 
eigenvectors,  p^ .  This  projection  generated  the  eigenvectors 
which,  in  a  least  squares  sense,  came  closest  to  matching  the 
desired  eigenvectors. 


In  this  manner,  an  aclevable  eigenvector  was  found 
corresponding  to  each  desired  eigenvalue,  and  the  P 
matrix,  consisting  of  columns  of  eigenvectors,  was  formed. 
Naturally,  having  to  use  projected  eigenvectors  instead  of 
the  desired  eigenvectors  caused  some  deviation  from  the 
optimal  handling  qualities  designed  for  in  Section  III.  This 
was  unavoidable  since  the  number  of  inputs  was  less  than  the 
number  of  states. 

Once  the  projected  eigenvectors  were  found,  the  MODES 
program  executed  a  straight-forward,  full-state  feedback 
e igenstructure  assignment  algorithm.  The  matrix 


A  -  PaAPa 


(4.3) 


has  already  been  found.  Each  column  of  P.  ,  considered  as  a 

A 

vector  p  ,  must  satisfy  the  basic  eigenvalue  problem 


A  P,.  =  \  Pl 


or  since  A  =  A  +  BK 


[A  +  BKJp^  =  V  p^ 

Rearranging  this  equation  leads  to 
[A  -  X  I  ]  p  ^  +  BKp  =  0 


or,  in  matrix  partitioned  form 

[A  -  X  I  IB]  f_P._l  =  0 

K 

This  is  tantamount  to  requiring  that 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


lie  in  the  null 


space  of  (A  -  X  I  |  B]  .  This  has  already  been  taken  care 
of  earlier  by  the  eigenvector  projection  process  and  the 
lower  partition  of  the  projected  eigenvector  determines  KP  . 
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It  remained  only  to  postmultiply  this  by  the  inverse  of  P4 
to  obtain  the  unique  feedback  gain  matrix  K. 


Application  of  MODES 

A  total  of  twelve  feedback  matrices  were  found  using  the 
MODES  program.  Three  sets  of  desired  eigenvalues  and 
eigenvectors  corresponding  to  good,  fair,  and  poor  handling 
qualities  were  used  at  each  of  the  four  flight  conditions. 

For  use  in  the  simulation  it  was  necessary  to  calculate 
the  feedback  gain  matrices  in  the  body  axis  frame.  The  bare 
airframe  X-29A  A  and  B  matrices  calculated  in  Section  II 
were  already  in  this  frame.  However,  the  optimal  handling 
quality  A  matrices  from  which  the  desired  e igenstructure 
were  obtained  were  found  in  the  stability  axis  frame  (this 
frame  was  necessary  since  most  of  the  handling  quality 
specifications  are  expressed  in  this  frame  [3]).  Therefore, 
it  was  necessary  to  transform  the  stability  axis  dimensional 
derivatives  found  in  Section  III  into  the  body  axis, 
reformulate  the  A  matrix,  and  decompose  the  new  optimal 
matrix  into  the  desired  eigenvalues  and  eigenvectors. 

Again,  a  computer  program  was  written  to  facilitate 
this.  The  program  allowed  the  user  to  enter  dimensional 
derivatives  in  either  of  the  reference  frames.  The  A  and 
B  matrices  were  then  formed  in  whichever  frame  was  desired. 
The  optimal  matrix  was  then  decomposed  into  its  eigenvalues 
and  eigenvectors  using  an  IMSL  routine. 

Using  MODES,  it  was  a  relatively  quick  process  to 
formulate  the  twelve  feedback  gain  matrices.  The  personal 
computer  version  of  MODES  was  used  on  both  an  IBM  PC  and  IBM 
AT.  For  the  Mach  =  0.6,  altitude  =  10k  ft  flight  condition, 
and  the  e igenstructure  corresponding  to  good  handling 
qualities,  the  feedback  gain  matrix  calculated  by  MODES  is 


shown  below: 


K 


2 . 5642E-03 
5. 4848E-04 
3 . 8463E-04 


-8 . 0207E+00 
2 . 6434E+00 
-4 . 0206E+00 


2 . 4117E+00 
-3 . 9755E-01 
2 . 2998E-01 


-7.1717E-03 
1 . 6116E-03 
-8 . 5830E-04 


As  an  aid  to  understanding  this  matrix,  the  following  example 
is  given: 

6  .  .  =  —  2 . 5642E-03  *  u  .  -  8.0207E+00  *  a 

s  t  r  a  k  e  p  *  r  t  pen 

+  2 . 4112E+00  *  q  .  -  7.1717E-03  *  6 

•pert  pert 

and  similarly  row  2  provides  the  flaperon  command  and  row  3 


provides  the  canard  command. 

MODES  also  formulated  the  A 

matrix: 

-1 . 09  26E-02 

6 . 7087E+01 

-3. 4515E+01 

-3. 2126E+01 

A  — 

-1. 2758E-04 

-1.3133E+00 

9.9998E-01 

-2 . 8178E-03 

A  - 

8 . 4952E-04 

-7.5379E+00 

-3 . 5836E+00 

5. 5559E-03 

0 . 0000E+00 

0 . 0000E+00 

1 . 0000E+00 

0 . OOOOE+OO 

Of  course,  the  critical  issue  was  whether  the  A  matrix 
still  possessed  the  dynamic  characteristics  designed  for  it. 

A  comparison  of  the  eigenvalues  of  the  A  matrix  to  those  of 
the  original  optimal  A  matrix  of  Section  III  revealed 

Achievable  A  Matrix 
Eigenvalues 

X  .  =  -2 . 447761E+00  ±  j ( 2 . 504523E+00 ) 

I  ..  i 

V  „  =  -6 . 168339E-03  ±  j ( 7 . 340013E-02 ) 

Optimal  A  Matrix 
Eigenvalues 

\t  ,  =  -2 . 447769E+00  ±  j (  2 . 504519E  +  00 ) 

X.  „  =  -6 . 212100E-03  ±  j  (  7 . 334608E-02 ) 
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As  can  be  seen,  a  very  good  match  was  made.  As  expected,  the 
eigenvectors  did  not  compare  numerically,  however  they  still 


may  contain  the  same  approximate  mode  shapes  for  which  they 

A* 

were  designed.  From  the  A  matrix,  an  approximate  value  was 
found  for  Z_,  =  848.14  .  This  was  used  to  calculate  values  for 

'JL 

n/ct  and  CAP  for  the  augmented  system: 

A  Optimal  A 

n/ct  26.36  26.39 

CAP  0.465  0.464 

Again,  a  very  good  match  was  made.  It  appears  as  though  the 
desired  e igenstructure  came  through  the  eigenstructure 
assignment  process  in  good  shape. 

Feedback  Matrix  Preparation 

In  order  for  the  pilot  to  be  able  to  maneuver  the 
aircraft  throughout  the  envelope,  the  elements  of  the 
feedback  gain  matrices  for  each  of  the  handling  quality 
levels  were  curve  fit  together  as  a  function  of  aircraft 
dynamic  pressure.  A  third-order  equation  was  found  to  be 
sufficient  and  yielded  an  error  of  less  than  2%  between  the 
exact  value  and  the  curve  fit  value.  The  twelve  equations 
were  evaluated  every  25  msec  in  the  digital  simulation  to 
provide  the  feedback  gains.  The  implementation  of  the 
modified  control  system  into  the  X-29A  simulation  will  be 
described  in  Section  VI.  First,  in  Section  V,  a  pilot-model 
analysis  will  be  described  which  was  used  to  predict  the 
handling  qualities  of  the  closed-loop,  augmented  system. 


V.  Pilot  Model  Analysis 


v 


■s? 


-> 


Background 

With  the  feedback  matrices  developed  in  the  last 
section,  the  complete  closed-loop  system  was  ready  to  be 
tested.  However,  before  the  new  flight  control  system  was 
evaluated  on  the  man-in-the-loop  simulator,  an  attempt  to 
predict  the  handling  qualities  of  the  aircraft  was  made.  To 
accomplish  this,  a  pilot-model  analysis  was  performed  and  the 
Neal-Smith  criteria  was  used  to  predict  pilot  handling 
quality  ratings  1161. 

In  the  early  1970's,  Neal  and  Smith  published  a  report 
which  investigated  the  effects  of  control  system  dynamics  on 
the  longitudinal  flying  qualities  of  tighter  aircraft.  They 
used  the  USAF  NT-33  variable  stability  aircraft  to  evaluate 
57  different  combinations  of  flight  condition,  control  system 
dynamics,  and  short-period  dynamics.  One  of  the  important 
products  of  this  study  was  a  pilot-model  analysis  tool  which 
enables  the  control  system  designer  to  predict  what  the  pilot 
handling  quality  ratings  will  be  [16:  Sec  VI]. 

To  apply  the  Neal-Smith  technique,  one  needs  the 
closed-loop  aircraft  system,  as  well  as  an  outer  loop  which 
the  pilot  will  close  (usually  a  theta  command  loop).  A 
typical  configuration  is  shown  in  Figure  7.  The  pilot  model 


was  briefly  discussed  in  Section  I.  Once  again,  Kp  is  the 

_  3  s 

pilot  gain,  or  how  much  he  moves  the  stick;  e  is  the 


time  delay  due  to  pilot  reaction  time,  neuromuscular  delay. 


and  other  effects;  x 


and 


p  ,, 


are  the  pilot  lead  and  lag 


time  constants,  respectively.  The  pilot  will  adjust  Kp  , 
t  ,  and  to  attain  certain  performance  standards. 


described  later  in  this  section 
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PILOT  MODEL 


AIRFRAME 


l  J  _  ,  Tp  s  +  * 

e  ^‘Kd  e  ,3s  — 1 - 

Tp,3  *  1 


Figure  7.  Block  Diagram  of  Pitch  Attitude  Tracking 

Nichol's  charts  are  used  extensively  in  the  pilot-model 
analysis  to  relate  the  open-loop  transfer  function,  B/8f,  to 
closed-loop  performance.  Using  the  model  in  Figure  7,  the 
following  design  criteria  were  proposed. 

Bandwidth  (BW)  .  Bandwidth,  by  the  Neal-Smith 
definition,  is  the  frequency  for  which  the  closed-loop  Bode 
phase  is  equal  to  -90  degrees.  It  is  a  measure  of  how 
quickly  the  pilot  can  move  the  airplane's  nose  toward  the 
target.  For  high-speed  flight,  Neal-Smith  recommend  a 
minimum  bandwidth,  of  3.5  rad/sec.  The  pilot  will 

adjust  his  gain  (Kp)  to  attain  this  bandwidth. 


Drool 


Droop  is  defined  as  the  maximum  excursion  of 


closed-loop  Bode  amplitude,  /8  I  ,  below  the  0  db  line 

c  ma  X 

for  frequencies  less  than  (BW  >min  •  For  a  reasonable  amount 
of  damping,  droop  is  a  measure  of  how  slowly  the  nose  settles 
down  on  the  target  (steady-state  tracking  error).  The  amount 


of  droop  which  the  pilot  will  allow,  according  to  Neal-Smith, 
Is  3  db  below  the  0  db  line. 


i 

P 

v“ 

C 

V 
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PIO  Tendency  .  The  tendancy  to  oscillate,  or  to 
generate  Pilot  Induced  Oscillations,  is  defined  in  terms  of 
the  Bode  magnitude  of  any  closed-loop  resonant  peak. 


/ 


9 


,  which  results  from  the  pilot  trying  to  meet  the 


performance  standards.  The  pilot  will  add  lead  to  decrease 
the  resonance  peak,  however  droop  will  increase.  As 
resonance  grows,  the  pilot  will  complain  first  of  overshoot, 
then  oscillations,  and  finally  PIO  for  high  resonance  values, 
|0/0  |  >  11  db. 

1  c  1  wax 

Pilot  Compensation  .  The  pilot's  workload  is  related  to 
the  amount  of  phase  lead  or  lag  he  is  adding  to  the  system. 

If  the  pilot  must  add  too  much  lead,  he  will  state  that  the 
aircraft  response  is  sluggish,  while  too  much  lag  generates 
comments  about  having  to  fly  the  aircraft  smoothly.  Pilots 
prefer  to  add  a  little  lead  rather  than  lag.  The  amount  of 
pilot  compensation  is  defined  by 


ARG 


pc 


ARG  ( 


+  1 


+  1 


W  =  (BW) 

'  rru-  n 


(5.1) 


In  summary,  one  needs  to  determine  the  values  of  Kp  , 
t  ,  and  x  which  will  minimize  l9/0  j  while 

P  l  P  i  C  TT.  a  X 

maintaining  a  minimum  bandwidth  of  3.5  rad/sec  and  a  maximum 

1  8 

droop  of  3  db.  The  values  of  ARG  and  I  /g  j  which 

result  can  be  used  to  predict  the  pilot's  handling  quality 
rating  as  shown  in  Figure  8.  This  figure  is  based  on  the 
experimental  results  of  the  Neal-Smith  study. 
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Figure  8.  Neal-Smith  Criterion  for  Fighter  Maneuvering  Dynamics 


The  methology  to  determine  K  ,  t  ,  and  z  will  be 

P  P  i  p  2 

shown  by  example  later  in  this  section.  First,  a  simplified 
criteria  developed  by  Neal-Smith  will  be  discussed  [16:  Sec 
8.2].  This  method  for  predicting  handling  qualities  does  not 
require  finding  the  values  for  K  ,  z  ,  and  x  .  The  method 

P  P  j  P  ; 

takes  advantage  of  the  fact  that  the  amount  of  pilot 

compensation  added,  ARGa(j  /  is  related  to  the  open-loop  phase 

of  the  uncompensated  pilot  plus  airplane  at  u  =  ^Bw^min  * 

Also,  the  slope  of  the  compensated  amplitude-phase  curve  in 

the  neighborhood  of  u  =  (BW)m^n  is  an  approximate  measure  of 

how  large  the  closed-loop  resonance  will  be.  Therefore,  a 

rough  guess  to  the  amount  of  pilot  compensation  required  and 

the  closed-loop  resonance  are  determined  by  the  parameters 
d  A 

ARGad  and  ( •  Equations  for  these  two  parameters 
are  presented  below  [16:  117,118]: 

ARGa<J  -  [  M>G(^-)jBWniin  -  17.2(BW)mln  ]  ,de9  (5.2) 


dA 

( - >ad 

dARG  aa 


[d  ||r-|  /  d  ( log  U))] 


(BV) 


min 


d(ARG-f-)  /  d ( log  w) 

5 


-  39 . 6 ( BW) 


(BW) 


min 


min 


(5.3) 


The  parameters  are  compared  on  Figure  9  to  obtain  the  pilot 
ratings.  One  of  the  objectives  of  this  part  of  the  study  was 
to  compare  the  predictions  of  the  handling  qualities  of  the 
same  system  by  using  the  simplified  criteria  and  the  full 
pilot  model  analysis.  It  will  be  seen  that  it  became 
necessary  to  use  the  simplified  criteria  for  certain 
high-order  systems. 


X-29A  FI Ight  control  System 

As  a  reminder,  the  primary  purpose  of  using  the 
Neal-Smith  criteria  in  this  study  was  to  predict  the  pilot 
ratings  which  would  be  given  when  evaluations  were  made  using 
the  man-ln-the-loop  simulation.  Therefore,  It  was  necessary 
lo  duplicate  the  control  system  which  would  be  used  on  the 
LAMARS  during  this  phase  of  the  study. 

For  the  LAMARS  simulation,  the  basic  design  used  in  the 
actual  X-29A  Analog  Reversion  (AR)  mode  would  be  kept,  the 
major  changes  being  the  deletion  of  the  proportional  plus 
integral  compensator  and  the  addition  of  full-state  feedback. 
For  the  Neal-Smith  pilot  analysis,  the  control  system  used  is 
shown  i>n  Figure  10.  This  will  be  referred  to  as  the 
high-order  system.  As  can  be  seen,  a  number  of  first-order 
lags,  as  well  as  fourth-order  actuator  models  were  Included. 

A  second-order  Pade '  approximation  was  used  to  model  the 

_  Oe 

pilot  time  delay,  e  .  These  higher-order  effects  were  not 
taken  into  account  when  the  feedback  matrices  were  developed 

and  will  likely  degrade  the  system  performance.  When  the 

□ 

high-order  system  open-loop  transfer  function  (  /0  )  was 

formed  using  TOTAL,  it  contained  a  23rd  order  numerator  and  a 
29th  order  denominator.  However,  what  became  a  more 
difficult  problem  was  that  the  coefficients  for  the 
polynomials  were  of  the  power  10iU-  10' 8  .  Due  to  the  large 
coefficients  of  the  high-order  system  and  the  fact  that  the 
exponential  notation  used  by  TOTAL  provided  only  4 
significant  digits  in  its  output  (i.e.  0.1234E+28),  accurate 
Nlchol's  charts  could  not  be  obtained  for  these  large 
transfer  functions  (specifically,  the  lower  frequency  part  of 
the  curve  was  lost).  For  this  reason,  it  was  decided  to 


apply  the  simplified  Neal-Smith  criteria  to  the  high-order 
'/p  transfer  function,  and  apply  the  full  Neal-Smith 


Figure  10.  X-29A  Higher-Order  Mathematical  Model 


analysis  to  the  low-order  system,  shown  in  Figure  11. 

The  Nichol's  charts  drawn  by  TOTAL  were  not  adequate  to 
apply  the  Neal-Smith  criteria  in  that  they  did  not  display 
the  constant  M  and  N  curves.  Therefore,  a  program  was 
written  which  drew  a  complete  Nichol's  chart  using  the 
DISSPLA  graphics  package  and  which  emphasized  the  Neal-Smith 
criteria  of  90  degrees  phase  for  bandwidth  and  3  db  for 
low-frequency  droop. 

Neal-Smith  Analysis 

Recall  that  feedback  matrices  were  formed  for  four 
flight  conditions  in  each  of  three  handling  quality  levels 
(Good,  Fair,  and  Poor).  The  flight  conditions  corresponding 
to  M  =  .6,  Alt  =  10K  for  each  level  (Cases  1,  5,  &  9)  were 
selected  to  be  evaluated  using  the  Neal-Smith  criteria. 

Case  1_.  The  8/5  transfer  functions  were  found  by  using 
the  resolvent  matrix  for  the  A  +  BK  system  and  were 
provided  by  the  MODES  program,  discussed  in  Section  IV.  It 
remained  only  to  use  TOTAL  to  combine  the  series  and  parallel 
paths  into  the  open-loop  transfer  functions  8/Fs  and  8/9e . 
The  Bode  plots  for  the  lower-order  (LO)  and  higher-order  (HO) 
systems  are  shown  in  Figures  12  and  13.  From  tabular 
printouts.  It  was  found  that  the  phase  margin  degraded  from 
95  for  lower-order  8/Fs,  to  58  when  the  pilot  pure  time 
delay  was  added  -  8/9f  (LO).  Similarly,  the  gain  margin  went 
from  15  db  to  4.25  db.  When  the  higher-order  effects  were 
taken  into  account,  the  phase  margin  for  9/Fr  dropped  to  65°, 
and  for  0/?e  to  18  .  Likewise,  the  gain  margin  was  7.5  db 
for  8/Fs  and  1.09  for  8/0e.  Clearly  the  stability  margins 
for  the  higher-order  transfer  function  were  decreasing. 
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gure  11.  X-29A  Lower-Order  Mathemat 


At  this  point  the  simplified  criteria  was  applied  to 

both  the  higher-order  and  lower-order  systems.  For  the  lower 

d  |S/F  ;  d  ( ARG  0/F  ) 

order  system,  aUftj-,T  -  -33  db/dec,  gfn)g-,.r-  =  -230^ 

and  (ARG  0/Fs)BW  =  -136  deg.  Therefore,  ARG  j  =  -196 
1  min 

degrees  and  (dA/dARG)a<^  =  +.0895  db/deg.  From  Figure  9,  this 
predicts  Level  2  handling  qualities.  Similarly  for  the 
higher-order  system,  ARGa(j  =  “200.2  deg  and 
(dA/dARG)a(j  =  +.0786  db/deg.  These  numbers  correspond  to 
Level  2/Level  3  handling  qualities  from  Figure  9.  Thus,  the 
simplified  criteria  Indicates  that  the  system  which  was 
designed  to  yield  good  handling  qualities  may  In  fact  be 
Level  2  or  worse.  Clearly  the  higher-order  effects  and  pilot 
time  delay  have  taken  their  toll  on  the  system  performance. 
Next,  ,  t  ,  and  t  were  found  for  the  lower-order 

system  using  the  method  described  Reference  16,  Section  6.6. 
First,  K  was  adjusted  to  bring  the  open-loop  amplitude-phase 

r 

curve  up  on  the  Nlchol’s  chart  so  that  it  crossed  the  0  db 
line  at  -180  degrees  of  phase.  Any  further  increase  in  Kp 
would  drive  the  system  unstable.  As  can  be  seen  from  Figure 
14,  a  value  of  1.63  accomplishes  this.  It  is  clear  that  lead 
compensation  by  the  pilot  can  be  used  to  bring  the  curve  to 
the  right  in  order  to  meet  the  bandwidth  criteria  (-90'  of 
phase  at  «  =  3.5  rad/sec).  As  pointed  out  by  Neal-Smith,  it 
is  desirable  to  have  the  phase  at  w  =  3.5  rad/sec  near  -130 
degrees  in  order  to  have  low  resonance  (16:55).  From  Figure 
14,  the  phase  at  w  =  3.5  rad/sec  was  -198' .  Therefore,  the 
pilot  would  have  to  provide  68"  of  phase  lead  at  *  -  3.5 
rad/sec.  From  Figure  15,  taken  from  the  Neal-Smith  report,  a 
phase  lead  of  68'  corresponds  to  Tf  ^  -  2.6  (xf  =  0),  or 

■  =  0.74  sec.  Therefore,  the  pilot  model  would  look  like 

1 
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Gp  =  ( 1 . 63 )e  ’ 3s ( 0 . 74s  +  1)  (5.4) 

However,  when  this  new  pilot  model  was  Included  in  the  9/9e 
transfer  function  and  plotted  on  the  Nichol's  chart,  the 
performance  standards  were  not  quite  met.  An  iterative 
process  began  until  the  correct  values  of  Kp  =  0.82  and 
x  =0.42  are  found. 

pi 

The  resulting  Nichol's  chart  and  closed-loop  Bode 

diagram  are  shown  in  Figures  16  and  17  respectively.  The 

amount  of  pilot  phase  added  at  w  =  3.5  rad/sec  was  56'  and 

|8/6r  |  =  5.4  db.  Therefore,  from  Figure  8,  the  aircraft 

would  exhibit  Level  2  handling  qualities,  as  the  simplified 

criteria  predicted.  This  gives  us  some  faith  in  appling  th‘e 

simplified  criteria  to  the  higher-order  systems. 

Case  jj .  Figures  18  and  19  show  the  Bode  plots  for  the 

lower  and  higher  order  systems,  respectively.  The  phase  and 

gain  margins  indicate  that  they  are  all  unstable  systems. 

Figure  20  shows  the  Nichol's  chart  for  the  lower-order  system 

with  a  value  for  Kp  of  0.35.  Several  attempts  were  made  to 

meet  the  performance  standards  using  pure  lead  (x  =0)  but 

P  2 

no  configuration  was  found.  When  both  x  and  x  were 

O  O 

r  i  ■  : 

adjusted,  slightly  better  results  were  obtained.  A  value  of 
K  =  0.06,  x  =  5.0,  and  x  =  0.005  was  used  to  produce  the 

P  P  ,  P  ; 

Nichol’s  chart  in  Figure  21  and  the  closed-loop  Bode  diagram 

in  Figure  22.  The  value  of  x  used  here  is  unreasonable. 

p. 

The  amount  of  pilot  phase  added  was  86”  and  19/9. !max  =  4.4 
db.  The  maximum  low-frequency  droop  was  6.25  db.  This  is 
probably  a  Level  3  system  as  seen  in  Figure  8. 

Appling  the  simplified  criteria  to  the  higher-order 
system  resulted  in  (ARG)a(j  =  -254  deg  and  (dA/dARG)acj  =  .18 
db/deg.  Figure  9  shows  this  to  be  Level  3,  almost  Level  2. 
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Figure  16.  Case  1,  Nichol's  Chart  With  Pilot  Equalization 
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Case  9.  Figures  23  and  24  are  the  Bode  plots  for  the 
lower  and  higher-order  systems,  respectively.  Again,  the 
phase  and  gain  margins  indicate  that  the  systems  are 
unstable.  Figure  25  is  the  Nichol's  chart  for  the 
lower-order  system  with  Kp  =  0.04.  It  is  obvious  that  a  good 
deal  of  compensation  Is  required  to  meet  the  performance 
standards.  A  pilot  model  with  K  =  0.03,  =  10.,  and 

t  =  0.001  gave  the  closest  agreement  with  the  performance 

P  2 

standards,  as  shown  In  Figure  26.  The  closed-loop  Bode 


diagram  is  shown  in  Figure  27.  The  amount  of  pilot  phase 
added  was  88  and  |0/8. Imax=  19.5  db.  The  maximum  droop  is 


2.5  db.  This  is  a  Level  3  system  with  strong  PIO  tendancles, 
as  seen  in  Figure  8. 

Appling  the  simplified  criteria  to  the  higher-order 
system  resulted  in  (ARG)a(j  =  -280  deg  and  (dA/dARG)ad  =  .23 
db/deg.  From  Figure  9,  this  is  a  Level  3  system. 


Summary 


The  addition  of  even  just  the  pilot  delay  and  stick 
prefilter  noticably  degrades  the  performance  of  the  systems. 
According  to  Neal-Smlth,  the  system  designed  to  be  Level  1  is 
Level  2;  the  system  designed  to  be  Level  2  is  almost  Level  3; 
and  the  system  designed  to  be  Level  3  is  very  Level  3, 
possibly  uncontrollable.  The  small  values  of  Kp  could  be 
helped  by  reducing  the  stick  gain  from  2.0  to  1.0.  It  was 
satisfying  to  see  good  correlation  between  the  simplified 
criteria  applied  to  the  higher-order  systems  and  the  full 
pilot  model  analysis,  applied  to  the  lower-order  systems. 

P 1  lot- in-the -loop  evaluations  will  be  examined  in  the  next 
section . 


*  f  ^  Wf  : 
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Figure  23.  Case  9,  Low-Order  Bode  Plots 
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VI.  Manned  Simulation  Analysis 


The  final  step  in  the  evaluation  process  was  to 
Integrate  the  full-state  feedback  control  system  into  the 
LAMARS  X-29A  simulation  and  conduct  man-ln-the-loop 
evaluations  for  a  pitch  tracking  task.  It  is  generally 
accepted  that  ground-based  man- in-the-loop  simulators  are  a 
very  useful  tool  but  should  never  be  considered  the  final 
word  on  the  acceptibility  of  a  control  system.  They  are  best 
utilized  in  the  role  of  bridging  the  gap  between  paper 
analysis  and  flight  test. 

The  LAMARS  simulation  of  the  X-29A  received  extensive 
scrutiny  from  government  agencies  and  industry  when  it  was 
used  as  the  Air  Force  safety-of-f light  evaluation  simulator. 
It  is  recognized  as  being  a  superior  simulation  and 
therefore,  the  results  of  the  current  study  should  be  valid 
within  the  restrictions  of  the  control  system  design.  The 
most  notable  restrictions  were  the  small  perturbation 
assumption  and  the  piecewise  linearization  of  feedback  gains 
throughout  a  limited  flight  envelope. 

Description  of  the  Flight  Control  Development  Laboratory 

The  Flight  Control  Development  Laboratory  ( FCDL )  is 
located  within  the  Flight  Control  Division  of  the  Flight 
Dynamics  Laboratory.  The  primary  mission  of  the  facility  is 
to  evaluate  flight  control  system  concepts  in  as  realistic  an 
environment  as  necessary  through  the  use  of  several 
man-ln-the-loop  simulators.  The  simulators  Include  three 
motion-base,  one  static-base,  and  one  manned  combat  station. 

A  brief  description  of  the  hardware  used  in  this  study  is 
presented  below  in  order  to  give  the  reader  a  sense  of  the 
level  of  fidelity  used  for  the  study. 
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Hardware .  The  LAMARS,  shown  in  Figure  28,  provides  high 
quality  flight  cues  in  a  five  degree-of -freedom  (lateral, 
vertical,  roll,  pitch,  and  yaw),  beam-type  motion  system. 

The  20-foot-diameter  dome  contains  a  single-seat  fighter 
cockpit  and  spherical  dome  display.  The  LAMARS  motion  system 
consists  of  a  30-foot-long  horizontal  beam,  gimbaled  at  the 
rear  and  driven  by  hydraulic  actuators  to  provide  ±  10  feet 
of  both  vertical  and  lateral  motion  to  the  cockpit.  The 
cockpit  gimbal  system  is  mounted  on  the  forward  end  of  the 
beam  and  provides  angular  rotation  of  t  25  degrees  in  pitch, 
roll,  and  yaw.  The  LAMARS  motion  system  performance  is 
summarized  in  Table  V. 

The  transient  motion  cues  provided  by  the  LAMARS  motion 
system  were  augmented  by  a  g-suit  which  is  programmed  to 
provide  positive,  sustained  cueing  for  load  factor  conditions 
above  one  gee. 

The  cockpit  was  configured  as  a  generic  fighter  cockpit 
with  a  standard  suite  of  up-front  instruments,  a  head-up 
display  (HUD),  throttle,  and  a  programmable  feel  system 
center-stick.  A  typical  LAMARS  cockpit  configuration  is 
shown  in  Figure  29.  On  the  inside  surface  of  the  dome  a  300 
degree  dynamic  sky-earth  image  was  projected  to  aid  the  test 
subject  in  visual  cues.  Also,  artificial  engine  noise 
provided  throttle  setting  feedback  to  the  test  subject  and 
added  to  the  overall  environment  realism. 

At  the  heart  of  the  simulation  are  four  Gould  S EL 
computers  which  communicate  to  each  other  through  1  MegaByte 
of  shared  memory.  The  computers'  CPU's  and  IPU's  are  run  in 
a  parallel  processor  fashion.  The  SEL's  are  32-bit  digital 
machines  with  high-speed,  real-time  and  scientific  computing 
capabilities.  A  real-time  hybrid  clock  provides  the  clock 
pulse  to  start  each  simulation  cycle  and  initiate  the 
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dlgital-to-analog  (DAC)  and  analog-to-dlgital  (ADC) 
conversions.  A  25  mllli-second  frame  rate  was  used  on  this 
simulation.  An  EAI  781  analog  computer  was  used  to  model  the 
Analog  Reversion  mode  of  the  X-29A  flight  control  system  and 
also  the  fourth-order  actuator  models.  A  CSPI  array 
processor  was  used  to  perform  table  look-ups  for  the 
aerodynamic  math  model.  Over  60,000  data  points  in 
breakpoint  format  were  used  in  the  data  package. 

An  Interactive  Machines  Inc.  (IMI)  455D  stroke  graphics 
system  was  used  to  generate  the  head-up  display,  shown  in 
Figure  30.  The  IMI  communicated  to  the  SELs  through  an 
Ethernet  system.  There  was  virtually  no  time  delay  in  the 
Ethernet/IMI  system  for  this  simple  HUD.  A  Megatek  Whizzard 
7000  calligraphic  system  was  used  to  provide  the  test 
conductor  with  a  real-time  display  of  pertinent 
simulator/alrcraf t  data  on  a  monitor  adjacent  to  the  test 
conductor's  station. 

The  Master  Control  Console  (MCC),  shown  in  Figure  31, 
provides  the  ability  for  one  simulation  operator  to  monitor 
and  control  the  entire  simulation.  A  complete  complement  of 
aircraft  Instruments  replicates  those  in  the  cockpit.  Video 
monitors  provide  Images  of  the  HUD,  LAMARS  cockpit 
over-the-shoulder,  and  LAMARS  bay.  A  joystick  and  set  of 
potentiometers  allows  the  simulation  operator  to  fly  the 
simulation  during  development  and  check-out.  Thirty  lighted, 
pushbutton  switches  provide  the  ability  to  control  various 
aircraft  subsystems,  to  switch  between  the  standard  AR  mode 
and  the  modified  control  system,  and  to  select  the  particular 
modified  system  to  evaluate  (Good,  Fair,  or  Poor  handling 
qualities).  Several  multichannel  Intercom  units  are 
available  to  provide  communication  between  the  test  subject, 
LAMARS  operator,  simulation  operator,  and  observers. 
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Software.  A  complete  description  of  the  software  used 


by  the  X-29A  simulation  is  beyond  the  scope  of  this  report. 


All  of  the  software  is  coded  in  FORTRAN  and  written  in-house. 


Relatively  few  changes  were  made  to  the  simulation  to  carry 
out  this  study  and  they  will  be  discussed  later  in  this 


section . 


Modification  of  the  X-29A  Longitudinal  Control  Path 


A  three-step  process  was  used  to  integrate  the  new 
control  system  into  the  simulation.  Before  discussing  this 
process,  an  explanation  needs  to  be  made  concerning  the  use 
of  the  word  "controller''.  In  this  report,  the  word 
controller  is  meant  to  describe  that  portion  of  the  flight 
control  system  which  interprets  the  stick  movements  and 
creates  four  state  perturbation  signals  for  the  twelve 
feedback  gains.  As  will  be  seen,  the  design  of  the 
controller  was  an  interesting  problem. 

The  pitch-rate  feedback,  proportional  plus  integral 
longitudinal  control  path  of  the  standard  X-29A  AR  mode  on 
the  analog  computer  was  replaced  by  the  a  simple  full-state 
feedback  system  with  no  controller  installed,  as  shown  in 
Figure  32.  Actually,  the  standard  AR  mode  remained  patched 
and  a  series  of  function  relays  activated  by  a  single  MCC 
button  was  used  to  switch  between  the  standard  AR  path  and 
the  modified  path. 

Initial  attempts  to  use  this  modified  control  path 
failed  to  keep  the  aircraft  from  diverging  in  pitch.  The 
cause  of  the  divergence  turned  out  to  be  the  added  delay 
introduced  by  the  fourth-order  actuator  models.  Recall  that 


the  actuator  time  constants  were  not  taken  into  account  when 


the  feedback  gains  were  calculated.  The  penalty  for  this  was 
made  apparent  by  the  diverging  aircraft.  When  the  actuator 


gure  32.  Modified  Analog  System,  No  Controller 


models  were  bypassed  on  the  analog  computer,  the  aircraft 
remained  very  stable.  One  of  the  primary  purposes  of  the 
flight  control  system  was  successfully  completed:  to  make  an 
unstable  aircraft  behave  in  a  stable  manner.  A  turbulence 
model  was  used  to  provide  five  ft/sec  RMS  of  simulated 
turbulence  and  the  aircraft  handled  it  quite  well,  as  seen  in 
the  time  histories  of  Figure  33.  This  showed  that  a  second 
advantage  of  feedback  had  been  realized:  reducing  the  effects 
of  disturbances . 

One  of  the  attractive  parts  of  the  modified  control 
system  was  that  essentially  only  three  summers  on  the  analog 
computer  were  used  to  accomplish  this  design  (keep  in  mind 
that  the  controller  was  not  built  yet).  The  feedback  gains, 
which  were  now  a  function  of  dynamic  pressure,  were 
calculated  on  the  digital  computer,  multiplied  by  the 
appropriate  state  perturbation  variable,  and  sent  out  to  the 
analog  on  DACs .  The  four  feedback  signals  for  each  of  the 
three  control  surfaces  were  summed  together  to  formulate  the 
control  surface  commands  (which  in  the  absence  of  the 
actuators  were  the  deflections  too).  These  three  control 
surface  commands  were  read  back  into  the  digital  computer 
through  ADCs.  The  commands  were  used  by  the  aerodynamic 
model  to  compute  the  aerodynamic  force  and  moment 
coefficients  for  the  equatlons-of-motion . 

In  view  of  the  simple  nature  of  the  aircraft 
stabilization  system,  it  was  decided  to  move  the  modified 
control  path  from  the  analog  computer  to  the  digital 
computer.  This  would  permit  a  more  flexible  system  to  be 
used  during  debugging  and  check-out,  allow  digital  recording 
of  control  system  variables,  and  alleviate  the  magnitude 
scaling  worries  of  the  analog  computer.  The  capability  to 
switch  between  the  standard  AR  mode  and  the  modified  path  was 


retained,  as  well  as  the  Interface  to  the  analog 

lateral/d irectlonal  control  path.  This  switch  to  the  digital 

computer  was  the  second  of  the  three-step  process . 

The  last  step  was  to  design  a  controller  which  would 
allow  a  pilot  to  maneuver  the  aircraft.  The  controller 
design  was  based  on  a  search  of  the  literature  and 
discussions  with  flight  control  engineers  in  the  Flight 
Control  Division.  Unfortunately,  most  work  being  done  in  the 
division  focuses  on  the  transfer  function  model,  for  Instance 
8 / S e .  The  design  of  the  controller  is  left  to  Industry. 
Likewise,  most  papers  on  control  system  design  seem  to  stop 
at  an  examination  of  the  transfer  functions.  For  this 
reason,  a  trial  and  error  approach  was  used  to  des.ign  the 
controller  . 

Experience  indicated  that  a  pitch-rate  command,  attitude 
hold  system  was  an  appropriate  controller  to  try.  Therefore, 
the  stick  was  made  to  provide  an  added  pitch-rate 
perturbation  for  that  particular  feedback  path.  However,  due 
to  the  large  feedback  gains  on  angle-of -attack,  this 
pitch-rate  commend  system  had  very  little  effect  on  the 
aircraft.  Clearly  the  alpha  feedback  was  providing  the 
static  stability  for  the  aircraft  and  it  would  need  some  form 
of  a  command  if  the  aircraft  was  to  be  successfully 
maneuvered . 

Recall  that  a2  =  w  -  U  q  or  az  =  U*(i  -  q).  This 

c  •?  c  g 

provided  an  equation  to  determine  an  appropriate  alpha 
command.  A  stick  gradient  for  both  normal  acceleration  and 
pitch-rate  were  devised  based  on  standard  X-29A  time 
histories  of  stick  pulses.  Thus,  a  stick  deflection  resulted 
in  both  a  normal  acceleration  command  (AZCOM)  and  pitch-rate 
command  (QCOM)  from  which  an  alpha-dot  command  was  calculated 


according  to  the  equation 

ALPDOTCOM  =  AZCOM/UFPS  +  QCOM  (6.1) 

ALPDOTCOM  was  integrated  to  give  an  alpha  command  to  be  used 
when  calculating  the  alpha  perturbation  signal.  The 
disadvantage  of  this  scheme  was  that  while  the  alpha-dot 
command  would  go  to  zero  when  the  stick  wasn't  deflected,  the 
alpha  command  would  remain  at  whatever  value  it  held  when  the 
rate  went  to  zero.  This  would  result  in  the  aircraft 
continuing  to  change  pitch  attitude  after  the  stick  was 
released.  This,  of  course,  was  undesirable.  Therefore, 
logic  was  introduced  into  the  controller  which  would 
integrate  the  alpha  command  to  zero  at  a  rate  of  2 
degrees-per-second  if  the  stick  was  not  deflected.  The  rate 
of  integration  was  found  by  experimentation  and  balanced 
attitude  overshoot  (small  rate)  with  pitch  bobble  (large 
rate).  Since  this  was  essentially  an  alpha  feedback  system, 
the  reference  for  orientation  stabilization  was  the  relative 
wind,  whereas  for  a  conventional  aircraft  it  would  be  the 
horizon.  This  resulted  in  a  control  system  which  handled 
disturbances  by  trying  to  maintain  a  constant  angle-of- 
attack . 

The  forward  velocity  and  theta  perturbation  command 
signals  were  less  of  a  problem  since  they  primarily  provided 
phugoid  stability.  In  fact,  with  the  strong  angle-of-attack 
feedback,  the  assumption  that  the  phugoid  motion  occurs  at 
constant  angle-of-attack  was  even  more  valid  than  usual.  It 
should  be  noted  that  forward  velocity  feedback  provided  very 
good  speed  stability.  However,  since  the  throttle  was  not 
one  of  the  controllers,  speed  stability  was  accomplished  by 
varying  the  total  drag  of  the  aircraft  through  control 
surface  deflections.  This  feature  became  unruly  when  it  came 
to  commanded  attitude  changes  and  the  capability  to  turn-off 


the  u  feedback  was  added  to  the  simulation. 

An  attitude  hold  feature  was  added  for  the  theta 
feedback.  This  was  accomplished  by  having  a  zero  theta 
perturbation  signal  while  the  stick  was  deflected;  when  the 
stick  was  centered,  theta  was  sampled  and  used  as  the 
reference  value  until  another  stick  deflection.  However, 
again  due  to  the  overwhelming  alpha  feedback,  attitude  hold 
was  only  marginally  successful. 

The  flight  control  system  appeared  as  shown  In  Figure 
34.  As  can  be  seen,  a  stick  prefilter  and  stick  gain  have 
been  added  so  that  the  system  more  closely  matches  that  used 
in  the  pilot-model  analysis.  Digitally  controlled  stick 
inputs  were  used  to  get  aircraft  time  histories.  The 
response  to  a  four  second,  one  inch  stick  input  is  shown  in 
Figures  35-37  for  each  of  the  feedback  gain  sets.  The  most 
notable  feature  is  the  slow  onset  of  pitch-rate  compared  to  a 
more  conventional  design.  This  is  because  the  system  is 
essentially  a  position  command  system  which  is  unusual  for 
up-and-away  flight.  The  difference  in  response  between  the 
three  designs  is  evident.  The  system  designed  to  yield  fair 
handling  qualities  is  more  sluggish  and  oscillatory  than  the 
system  designed  to  yield  good  handling  qualities.  As  can  be 
seen,  the  control  system  designed  to  give  poor  handling 
qualities  is  divergent  in  phugoid,  as  its  eigenvalues 
predicted . 

Test  Procedure 

In  order  to  get  the  best  correlation  possible  between 
the  pilot-model  analysis  and  the  pilot-ln-the-loop  analysis, 
the  same  pitch  tracking  task  used  in  the  Neal-Smith  study  was 
used  in  this  simulation  [16:231.  Therefore,  a  series  of 
pitch  pulses  lasting  110  seconds  was  programmed  and  used  to 


Figure  34.  Final  Flight  Control  System  Configuration 
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Figure  37.  "Poor"  System  Stick  Pulse  Response  (1  of  2) 
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drive  the  pitch  command  bar  on  the  HUD.  It  was  the  test 
subjects'  task  to  zero  out  the  error  between  the  command  bar 
and  the  pipper.  In  order  to  help  the  test  subject  remain  in 
the  longitudinal  axis  only/  the  lateral/directional 
equat ions-o£-motion  were  fixed.  In  addition,  the  test 
subject  was  told  to  leave  the  throttle  in  the  trim  position 
throughout  the  experiment.  As  a  result,  the  test  subjects 
could  concentrate  on  the  pitch  tracking  task  alone. 

A  test  matrix  was  formed  which  randomized  the  two 
control  systems  to  be  evaluated  (the  third  was  divergent). 

As  a  control  for  the  experiment,  all  the  test  subjects 
evaluated  the  standard  AR  mode.  This  provided  a  baseline  by 
which  to  judge  the  test  subjects.  Next  the  subjects  flew  a 
modified  system  with  the  u  feedback  engaged  to  show  that  this 
configuration  was  not  appropriate.  Then  the  subjects  flew 
each  of  the  modified  systems  twice  and  also  each  modified 
system  in  turbulence,  all  without  the  u  feedback.  The  test 
matrix  used  is  shown  in  Table  VI.  The  test  subjects  were  not 
told  which  system  they  were  evaluating. 


Table  VI 
Test  Matrix 


Subject  1 

Subject  2 

Subject  3 

Subject  4 

Subject  5 

AR 

AR 

AR 

AR 

AR 

* 

_  _  * 

* 

* 

* 

GOOD 

FAIR 

GOOD 

FAIR 

GOOD 

GOOD 

FAIR 

GOOD 

FAIR 

GOOD 

FAIR 

GOOD 

FAIR 

GOOD 

GOOD 

FAIR 

GOOD 

GOOD 

GOOD 

FAIR 

GOOD 

FAIR 

FAIR 

FAIR 

FAIR 

** 

** 

*  * 

*  * 

** 

FAIR 

FAIR 

FAIR 

GOOD 

GOOD 

** 

** 

** 

** 

** 

GOOD 

GOOD 

GOOD 

FAIR 

FAIR 

*  -  flown 

with  speed 

stability  engaged 

**  -  flown 

in  moderate 

turbulence 

After  each  test  point,  the  test  subject  was  asked  to 
give  a  Cooper-Harper  rating  (Figure  38),  a  PIO  rating  (Figure 
39),  and  any  other  comments.  The  comments  were  recorded  by  a 
tape  deck  Interfaced  into  the  Intercom  system  so  that  they 
could  be  analyzed  at  a  later  time.  Also,  digital  recording 
of  pertinent  simulation  variables  was  made  for  each  run.  The 
results  of  the  simulation  study  will  be  presented  in  the  last 
section . 


<*  A 


DESCRIPTION 

NUMERICAL 

RATING 

NO  TENDENCY  FOR  PILOT  TO  INDUCE  UNDESIRABLE 

MOTIONS 

1 

UNDESIRABLE  MOTIONS  TEND  TO  OCCUR  WHEN 

PILOT  INITIATES  ABRUPT  MANEUVERS  OR  ATTEMPTS 

TIGHT  CONTROL.  THESE  MOTIONS  CAN  BE  PREVENTED 

OR  ELIMINATED  BY  PILOT  TECHNIQUE. 

2 

UNDESIRABLE  MOTIONS  EASILY  INDUCED  WHEN  PILOT 

INITIATES  ABRUPT  MANEUVERS  OR  ATTEMPTS  TIGHT 

CONTROL.  THESE  MOTIONS  CAN  BE  PREVENTED  OR 

ELIMINATED  BUT  ONLY  AT  SACRIFICE  TO  TASK  PER¬ 
FORMANCE  OR  THROUGH  CONSIDERABLE  PILOT 

ATTENTION  AND  EFFORT. 

3 

OSCILLATIONS  TEND  TO  DEVELOP  WHEN  PILOT  INITIATES 

ABRUPT  MANEUVERS  OR  ATTEMPTS  TIGHT  CONTROL  . 

PILOT  MUST  REDUCE  GAIN  OR  ABANDON  TASK  TO 

RECOVER. 

4 

DIVERGENT  OSCILLATIONS  TEND  TO  DEVELOP  WHEN 

PILOT  INITIATES  ABRUPT  MANEUVERS  OR  ATTEMPTS 

TIGHT  CONTROL.  PILOTMUST  OPEN  LOOP  BY  RELEASING 

OR  FREEZING  THE  STICK. 

5 

DISTURBANCE  OR  NORMAL  PILOT  CONTROL  MAY ' 

CAUSE  DIVERGENT  OSCILLATION  .PILOT  MUST  OPEN 

CONTROL  LOOP  BY  RELEASING  OR  FREEZING  THE 

STICK. 

6 

Figure  39. 


Pilot  Induced  Oscillation 
(PIO)  Rating  Scale 
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VII.  Results  and  Recommendations 


LAMARS  Results 

The  man- in-the-loop  testing  on  the  LAMARS  was  the  final 
evaluation  tool  used  In  this  study.  Five  test  subjects  were 
selected  with  flying  experience  ranging  from  an  Air  Force 
Test  Pilot  School  graduate  to  an  aerobatics  private  pilot 
working  in  the  flight  controls  field.  Four  of  the  test 
subjects  had  previous  experience  flying  the  LAMARS  and 
evaluating  systems  on  it,  and  thus  were  more  accustomed  to 
the  sensory  feedbacks  the  LAMARS  provides.  All  testing  was 
accomplished  with  the  pilot  ’’blind"  to  the  particular  control 
system  being  evaluated.  They  were  also  unaware  of  total 
number  of  different  control  systems  being  evaluated. 

An  average  Cooper-Harper  rating  and  PIO  rating  is  given 
in  Table  VII  for  the  standard  AR  mode,  the  control  system 
designed  to  yield  good  handling  qualities,  the  control  system 
designed  to  yield  fair  handling  qualities,  and  the  runs  with 
turbulence . 


Table  VII 

LAMARS  Cooper-Harper  and  PIO  Ratings 


Standard 

AR 

Good 

System 

Fair 

System 

Turbi 

Good 

ilence 

Fair 

Cooper- 

Harper 

2.6 

7 . 5 

4 . 3 

7 . 5 

4.0 

Standard 

Deviation 

2.56 

1.18 

0.82 

-  -■  -  .  -  —  - 

1.0 

1.73 

PIO 

1.4 

3  .  4 

1.9 

3.7 

2.0 

Standard 

Deviation 

0.55 

0.73 

0.88 

0.96 

0.0 

Averaging  the  Cooper-Harper  ratings  Is  valid  as  long  as 
the  ratings  are  all  given  In  the  same  handling  quality  level, 
which  was  the  case  for  this  study.  The  ratings  for  each 
configuration  were  very  consistent,  as  seen  by  the  standard 
deviation.  Those  pilots  who  flew  the  same  conf Iguration 
back-to-back  showed  a  definite  learning  curve  and  usually 
gave  a  better  rating  the  second  time.  The  ratings  for  the 
standard  X-29A  AR  mode,  which  reflects  years  of  development 
put  into  the  control  system,  provided  a  baseline  by  which  to 
compare  the  modified  control  systems. 

Comments  on  the  "Good"  System.  As  an  overall  system, 
this  control  law  was  rated  as  Level  III  with  significant  PIO 
tendancies.  However,  a  close  review  of  the  pilot  comments 
recorded  during  the  testing  revealed  that  for  fine  tracking 
the  pilots  liked  the  system.  Only  during  the  second  half  of 
each  test  run,  when  gross  maneuvers  were  required,  did  the 
complaints  of  PIO  surface.  Some  typical  comments  were: 

"precise  tracking  is  better... I  get  into  a  PIO  on  gross 
maneuvers  which  are  almost  unstable" 

"large  Inputs  are  bad... have  to  add  lead  to  get  adequate 
performance" 

"sluggish  getting  star ted ...  large  captures  are  PIO  prone" 

When  asked  to  give  a  separate  rating  for  the  fine 
tracking  and  the  gross  acqusitlons,  the  fine  tracking  was 
given  a  Cooper-Harper  rating  in  the  3-4  range,  whereas  the 
gross  acquisition  was  given  a  rating  in  the  7-9  range,  with 
the  overall  rating  being  the  worse  of  the  two.  This 
disparity  between  fine  tracking  and  gross  maneuvering  was 
discussed  recently  In  a  paper  describing  the  AFTI/F-16 
control  modes  [21].  In  this  paper,  it  was  shown  that  a 


deadbeat  pitch  rate  response  is  best  for  fine  tracking  but 
lacks  the  quick  normal  acceleration  response  needed  for  gross 
acquisition.  The  pitch-rate  response  for  the  modified 
control  system  allowed  very  little  pitch-rate  overshoot  and 
in  addition,  the  normal  acceleration  build-up  was  slow.  This 
is  believed  to  have  contributed  to  the  PIO  problem  during 
large  Inputs. 

Another  possible  source  for  the  PIO  was  that  the  stick 
was  more  sensitive  than  the  other  control  system  evaluated. 

In  other  words,  the  system  was  running  with  a  high  gain.  The 
value  for  the  stick  gain  was  selected  as  a  compromise  between 
the  best  values  for  the  good  and  fair  systems.  It  was  felt 
that  keeping  the  stick  gain  constant  between  the  two 
configurations  was  necessary  for  comparison.  However,  this 
may  have  been  too  restrictive  and  a  smaller  stick  gain  should 
have  been  used  on  the  good  system. 

The  portion  of  the  control  system  to  blame  for  these 
short-comings  is  the  controller,  not  the  feedback  gains 
arrived  at  by  e lgenstructure  assignment.  As  mentioned  in 
Section  VI,  the  controller  design  was  not  conventional  for  an 
air-to-air  tracking  task.  It  is  felt  that  the  onset  of 
normal  acceleration  was  slow  contributing  to  the  PIO 
problems.  One  possible  solution  to  the  controller  problem 
would  be  to  force  the  feedback  for  static  stability  to  be  on 
pitch-rate  Instead  of  angle-of -attack  by  using  output 
feedback.  Then  a  controller  could  easily  be  devised  to 
create  a  pitch-rate  command  which  would  be  sufficient  to 
maneuver  the  aircraft. 

As  seen  by  the  Cooper-Harper  ratings,  the  addition  of  a 
moderate  level  of  turbulence  did  not  significantly  degrade 
the  performance  of  the  pilot.  This  was  to  be  expected  since 
state  feedback  should  damp-out  unwanted  disturbances. 
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Comments  on  the  "Fair"  System.  This  system  was  rated 
Level  II,  which  matched  the  handling  quality  level  for  which 
it  was  designed.  Some  of  the  comments  made  during  the 
testing  by  the  pilots  were: 

"oscillatory. . .precision  control  difficult" 

"more  controllable  for  large  commands" 

"annoying  bobbles  when  I  stop  the  stick" 

As  expected,  some  of  the  comments  centered  on  the  poor 
damping  of  the  overall  system. 1  However,  better  control  over 
gross  acquisition  kept  the  Cooper -Harper  ratings  in  the 
correct  level  for  a  system  having  minor  deflclences.  It  is 
believed  that  the  reason  this  system  did  not  experience  the 
PI 0  problems  of  the  good  system  was  that  the  stick  was  less 
sensitive.  This  resulted  in  a  response  which  was  more 
predictable.  Again,  turbulence  did  little  to  effect  task 
performance  for  the  fair  control  system. 

As  mentioned  in  an  earlier  section,  the  system  designed 
to  give  poor  handling  qualities  was  nearly  divergent,  even 
without  the  actuator  models.  Therefore  this  control  system 
was  not  evaluated  by  the  test  subjects. 

Conclusions 

The  culmination  of  this  thesis  project  can  be  summarized 
in  four  key  results: 

1)  Optimal  A  matrices  can  be  found  for  the  longitudinal 
axis.  This  matrix  can  be  weighted  toward  specific  tasks 
or  handling  quality  levels  by  the  judicious  selection  of 
the  handling  quality  parameters. 
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2)  Eigenstructure  assignment  with  four  states  and  three 
control  surfaces  produced  an  A  matrix  which  had  nearly 
identical  dynamic  character  1st ics  as  the  desired  A 
matrix.  The  necessary  projection  of  the  desired 
eigenvectors  onto  the  range  space  of  the  A  matrix  did 
not  change  the  eigenvectors  dramatically. 

3)  The  Neal-Smith  pilot-model  analysis  revealed  that  pilot 
delays,  actuator  models  (especially  fourth-order),  and 
control  system  dynamics  will  degrade  the  overall  system 
performance  by  approximately  one  handling  quality  level 
if  these  effects  are  not  taken  into  account  when, 
calculating  the  feedback  gains. 

4)  The  LAMARS  simulation  backed-up  the  Neal-Smith 
prediction  and  the  actuator  models  were  removed  to  get 
the  system  performance  back  closer  to  that  for  which  it 
was  designed.  The  type  of  controller  used,  in  this 
case  an  attitude  command/attitude  hold  system,  had  a 
significant  effect  on  the  handling  quality  ratings 

For  fine  tracking,  the  Cooper-Harper  ratings 
closely  matched  those  for  which  the  system  was  designed 
and  moderate  levels  of  turbulence  had  no  significant 
effect . 

Keep  in  mind  that  no  compensation  other  than  full-state 
feedback  through  a  simple  gain  was  used  to  attain  these 
results.  It  is  felt  that  the  technique  shows  great  promise. 

Recommendations 

A  number  of  areas  examined  during  this  project  warrant 
a  more  comprehensive  study.  In  particular,  the 


man- ln-the-loop  analysis  was  abbreviated  due  to  simulator 
availability.  A  careful  analysis  of  the  simulation  could 
easily  be  the  basis  for  a  follow-on  thesis.  It  would  be 
Interesting  to  apply  some  system  identification  techniques  to 
the  simulation  to  determine  exactly  what  the  total  system 
dynamics  are.  A  lower-order  equivalent  system  approach  could 
be  used  by  examining  some  LAMARS  time  histories  and 
determining  whether  the  system  displayed  conventional 
dynamics.  A  frequency  analyzer  could  be  used  to  obtain  Bode 
plots  of  the  total  system.  This  would  show  the  system  gain, 
bandwidth,  phase  margin,  and  gain  margin. 

A  different  type  of  controller  .could  be  developed  which 
may  yield  better  results.  In  particular,  the  stick  gain 
could  be  optimized  for  specific  tasks,  or  a  non-linear  stick 
gearing  used  to  get  different  response  rates.  Or,  as  in  the 
AFTI/F-16  paper,  a  blended  controls  mechanism  could  be 
employed  which  would  require  a  separate  prefilter  for  the 
normal  acceleration  and  alpha  command  paths.  Each  prefilter 
would  have  its  own  gain  which  would  vary. 

Also,  as  mentioned  earlier,  output  feedback  could  be 
used  to  force  the  pitch-rate  feedback  to  contain  the  static 
stability  augmentation.  This  would  allow  a  more  conventional 
controller  to  be  developed  which  would  permit  more  control 
over  pitch-rate  and  normal  accelearation  commands  and 
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The  use  of  eigenstructure  assignment  techniques  has  received  wide 
attention  as  a  tool  for  designing  flight  control  systems  for  aircraft 
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Therefore.,  development  of  a  method  for  choosing  the  desired 
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As  an  example-,  this  process  was  performed  on  a  model  of  the  X-29A 
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